{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "PoreP.ipynb",
      "provenance": [],
      "authorship_tag": "ABX9TyM2g8LxZUSzzqHDuH48JTIe",
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/Divyanshu-ISM/Oil-and-Gas-data-analysis/blob/master/PoreP.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "B7bBqlEX6nFJ",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "import numpy as np\n",
        "import pandas as pd\n",
        "import matplotlib.pyplot as plt\n",
        "%matplotlib inline"
      ],
      "execution_count": 1,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "ImyB3NIl6y0L",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "d = np.arange(20000,-1000,-1000)\n",
        "\n",
        "rhob = np.array([2.48,2.47,2.46,2.45,2.44,2.43,2.41,2.40,2.38,2.37,2.35,2.33,2.29,2.27,2.24,2.19,2.16,2.11,2.06,2.02,1.95])\n",
        "\n",
        "rhog = 2.65 #g/cm3\n",
        "\n",
        "rhof = 1.1 \n",
        "\n",
        "phi = (rhog - rhob)/(rhog - rhof)\n"
      ],
      "execution_count": 11,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "-YMtNTV_7ho6",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 545
        },
        "outputId": "d669841d-07ca-4c23-9f75-c353461d391c"
      },
      "source": [
        "plt.style.use('default')\n",
        "plt.figure(figsize=(3,6))\n",
        "\n",
        "plt.plot(rhob,d,label='RHOB')\n",
        "plt.plot(phi,d,label='PHI')\n",
        "\n",
        "plt.ylim(20000,0)\n",
        "# plt.xlabel('bulk density')\n",
        "plt.ylabel('Depth')\n",
        "plt.legend()\n",
        "\n",
        "plt.title('Porosity Determination')\n",
        "plt.grid()"
      ],
      "execution_count": 15,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUUAAAIQCAYAAAAW33yEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXwTZf4H8E/u9L5vSymtFFoKCkitcmppgcouHriCyiGH7oLKIWJdUY51QVhAVkDwp1JXQQVX0QUWDGdhWxQLpVsU5GhhQVKO3leSJs/vj3SmDb3TpjNJvu/Xq68kM09mnpmWL889EsYYAyGEEACAVOgMEEKImFBQJISQBigoEkJIAxQUCSGkAQqKhBDSAAVFQghpgIIiIYQ0QEGREEIaoKBICCENUFAkNpGeng6JRIKCggKhs+KwhL7HU6ZMQffu3QU5ty1RULQD3B8/96NWq9GzZ0/Mnj0bhYWFQmevzTZu3Ij09PROPWZBQYHFvVEoFPD398cDDzyA119/HVeuXLH62FVVVVi8eDEOHz7ceRm2M7/99hsWL16MnJwcobPSZSQ091n80tPTMXXqVCxduhSRkZGoqanBsWPH8OmnnyIiIgJ5eXlwdXUVOpsWjEYjDAYDVCoVJBIJAKBPnz7w9/fv1CBTUFCAyMhITJgwAWPGjIHJZEJxcTFOnDiBr7/+GhKJBB999BGeeuqpdh/71q1bCAgIwFtvvYXFixd3Wp47S1P3uLP99NNPuO+++7BlyxZMmTLFYp/BYIDJZIJKpbLJuYUiFzoDpO1Gjx6NgQMHAgCmT58OPz8/rFmzBt9++y0mTJhg9XEZY6ipqYGLi0tnZRUymQwymazTjtea/v3745lnnrHYdvnyZSQnJ2Py5Mno3bs3+vXr12X5aUllZSXc3Nw6fJyuvsd3UigUgp3blqj6bMceeughAEB+fj4AoLa2FsuWLUNUVBRUKhW6d++O119/HTqdzuJ73bt3xyOPPIJ9+/Zh4MCBcHFxwebNmwEAly5dwvjx4+Hr6wtXV1fcf//92L17d6Nzv/fee4iLi4Orqyt8fHwwcOBAbNu2jd9/Z3tX9+7dcebMGRw5coSv6g4fPhyXLl2CRCLB2rVrG50jMzMTEokEn3/+uVX3JyIiAunp6dDr9Vi5cqXFvpKSEsyZMwfh4eFQqVSIjo7GO++8A5PJBMBcAg0ICAAALFmyhM9zwxLj2bNn8cQTT8DX1xdqtRoDBw7Ed999Z3Ee7j4cOXIEf/rTnxAYGIi77roLADB8+HD06dMHubm5GDZsGFxdXREdHY2vvvoKAHDkyBEkJCTAxcUFMTEx2L9/f5PHbtimyP1ujx07hkGDBkGtVqNHjx74xz/+YfHdoqIivPLKK4iPj4e7uzs8PT0xevRonD59mk9z+PBh3HfffQCAqVOn8veAawJpqk2xsrIS8+fP5+9rTEwM/va3v+HOCqlEIsHs2bOxc+dO9OnTByqVCnFxcdi7d2+Tv8suxYjobdmyhQFgJ06csNi+bt06BoBt2rSJMcbY5MmTGQD2xBNPsA0bNrBJkyYxAGzcuHEW34uIiGDR0dHMx8eHvfbaa2zTpk3s0KFDTKvVsqCgIObh4cH+/Oc/szVr1rB+/foxqVTKvv76a/77H3zwAX+ezZs3s3Xr1rFp06axl156qVGe8/PzGWOMffPNN+yuu+5ivXr1Yp9++in79NNP2ffff88YY+zBBx9kAwYMaHTdf/rTn5iHhwerrKxs9t7k5+czAGzVqlXNpomKimIBAQH858rKSta3b1/m5+fHXn/9dbZp0yY2adIkJpFI2Msvv8wYY6yiooK9//77DAB79NFH+TyfPn2aMcZYXl4e8/LyYrGxseydd95h69evZ0OHDmUSicTiXnH3ITY2lg0bNoy99957bMWKFYwxxoYNG8ZCQ0NZeHg4W7BgAXvvvfdYbGwsk8lk7IsvvmDBwcFs8eLF7N1332VhYWHMy8uLlZWVNXuPud9tTEwMCwoKYq+//jpbv34969+/P5NIJCwvL49Pd+LECRYVFcVee+01tnnzZrZ06VL+HNeuXWOMMabVatnSpUsZADZz5kz+Hly8eJExZv57i4iI4I9pMpnYQw89xCQSCZs+fTpbv349Gzt2LAPA5syZY/E7AcD69evHQkJC2LJly9i7777LevTowVxdXdmtW7ea/V12BQqKdoD749+/fz+7efMm+9///se++OIL5ufnx1xcXNjVq1dZTk4OA8CmT59u8d1XXnmFAWAHDx7kt0VERDAAbO/evRZp58yZwwCwo0eP8tvKy8tZZGQk6969OzMajYwxxn7/+9+zuLi4NuW54T/YuLg4NmzYsEZpN2/ezACwX375hd+m1+uZv78/mzx5covnaUtQ/P3vf88AsNLSUsYYY8uWLWNubm7s119/tUj32muvMZlMxq5cucIYY+zmzZsMAHvrrbcaHfPhhx9m8fHxrKamht9mMpnYAw88wO6++25+G3cfBg8ezGpray2OMWzYMAaAbdu2jd929uxZBoBJpVJ2/Phxfvu+ffsYALZly5ZGx74zKAJgGRkZ/LYbN24wlUrF5s+fz2+rqanhf5+c/Px8plKp2NKlS/ltJ06caHRezp1BcefOnQwA+8tf/mKR7oknnmASiYRduHCB3waAKZVKi22nT59mANh7773X6FxdiarPdiQpKQkBAQEIDw/HU089BXd3d3zzzTcICwvDnj17AADz5s2z+M78+fMBoFEVODIyEikpKRbb9uzZg0GDBmHw4MH8Nnd3d8ycORMFBQX4+eefAQDe3t64evUqTpw40SnX9eSTT0KtVmPr1q38tn379uHWrVuN2gmt4e7uDgAoLy8HAOzYsQNDhgyBj48Pbt26xf8kJSXBaDQiIyOjxeMVFRXh4MGDePLJJ1FeXs5///bt20hJScH58+dx7do1i+/MmDGjyfY/d3d3i06gmJgYeHt7o3fv3khISOC3c+8vXbrU6vXGxsZiyJAh/OeAgADExMRYfFelUkEqNf/zNxqNuH37Ntzd3RETE4OTJ0+2eo6m7NmzBzKZDC+99JLF9vnz54Mxhn//+98W25OSkhAVFcV/7tu3Lzw9Pdt0jbZEHS12ZMOGDejZsyfkcjmCgoIQExPD/2FfvnwZUqkU0dHRFt8JDg6Gt7c3Ll++bLE9MjKy0fEvX75s8Q+R07t3b35/nz59sHDhQuzfvx+DBg1CdHQ0kpOTMXHiRDz44INWXZe3tzfGjh2Lbdu2YdmyZQCArVu3IiwsjG837YiKigoAgIeHBwDg/PnzyM3N5dsM73Tjxo0Wj3fhwgUwxrBo0SIsWrSo2WOEhYXxn5u63wBw1113Neo59vLyQnh4eKNtAFBcXNxi3gCgW7dujbb5+PhYfNdkMmHdunXYuHEj8vPzYTQa+X1+fn6tnqMply9fRmhoKH+fOQ3/ftqbTyFQULQjgwYN4nufm9PWoRkd6Wnu3bs3zp07h127dmHv3r345z//iY0bN+LNN9/EkiVLrDrmpEmTsGPHDmRmZiI+Ph7fffcd/vSnP/FBvyPy8vIQGBgIT09PAOaAMHLkSLz66qtNpu/Zs2eLx+M6Y1555ZVGpW3Onf85NXe/m+s9bm47a8MIurZ8969//SsWLVqE5557DsuWLYOvry+kUinmzJnDX5+tdeQabYmCooOIiIiAyWTC+fPn+f+ZAaCwsBAlJSWIiIho0zHOnTvXaPvZs2f5/Rw3Nzf84Q9/wB/+8Afo9Xo89thjePvtt5GWlga1Wt3k8VsK2KNGjUJAQAC2bt2KhIQEVFVV4dlnn201z63JysrCxYsXLarhUVFRqKioQFJSUovfbS6/PXr0AGAektLaMcTqq6++wogRI/DRRx9ZbC8pKYG/vz//uT3jHyMiIrB//36Ul5dblBab+vsRM2pTdBBjxowBALz77rsW29esWQMASE1NbdMxfvzxR2RlZfHbKisr8cEHH6B79+6IjY0FANy+fdvie0qlErGxsWCMwWAwNHt8Nzc3lJSUNLlPLpdjwoQJ2L59O9LT0xEfH4++ffu2mueWXL58GVOmTIFSqcSCBQv47U8++SSysrKwb9++Rt8pKSlBbW0tAPAD4u/Mc2BgIIYPH47Nmzfj+vXrjY5x8+bNDuW7K8hkskYlsh07djRqC+XGUzb3e2tozJgxMBqNWL9+vcX2tWvXQiKRYPTo0R3MddegkqKD6NevHyZPnowPPvgAJSUlGDZsGH788Ud88sknGDduHEaMGNHqMV577TV8/vnnGD16NF566SX4+vrik08+QX5+Pv75z3/yVdnk5GQEBwfjwQcfRFBQEH755ResX78eqampjdqTGhowYADef/99/OUvf0F0dDQCAwMt2gwnTZqEv//97zh06BDeeeeddl3/yZMn8dlnn8FkMqGkpAQnTpzAP//5T0gkEnz66acWAXbBggX47rvv8Mgjj2DKlCkYMGAAKisr8d///hdfffUVCgoK4O/vDxcXF8TGxuLLL79Ez5494evriz59+qBPnz7YsGEDBg8ejPj4eMyYMQM9evRAYWEhsrKycPXqVYvxfmL0yCOPYOnSpZg6dSoeeOAB/Pe//8XWrVv5UjAnKioK3t7e2LRpEzw8PODm5oaEhIQm20jHjh2LESNG4M9//jMKCgrQr18/fP/99/j2228xZ84ci04VUROw55u0UXPjFO9kMBjYkiVLWGRkJFMoFCw8PJylpaVZDBthzDxsIzU1tcljXLx4kT3xxBPM29ubqdVqNmjQILZr1y6LNJs3b2ZDhw5lfn5+TKVSsaioKLZgwQJ+yEvDPDccLqLVallqairz8PBgAJocnhMXF8ekUim7evVqK3fFjBuSw/3I5XLm6+vLEhISWFpaGrt8+XKT3ysvL2dpaWksOjqaKZVK5u/vzx544AH2t7/9jen1ej5dZmYmGzBgAFMqlY2G51y8eJFNmjSJBQcHM4VCwcLCwtgjjzzCvvrqq0b3oanf3bBhw5oc2tTc7wcAmzVrVqNj3zkkp6nvDhs2zOJ+19TUsPnz57OQkBDm4uLCHnzwQZaVldUoHWOMffvttyw2NpbJ5XKL4Tl3DslhzHxf586dy0JDQ5lCoWB33303W7VqFTOZTC1eS8P8tzYMy9Zo7jMRlXvvvRe+vr44cOCA0FkhToraFIlo/PTTT8jJycGkSZOEzgpxYlRSJILLy8tDdnY2Vq9ejVu3buHSpUvN9mATYmtUUiSC++qrrzB16lQYDAZ8/vnnFBCJoKikeIcNGzZg1apV0Gq16NevH9577z0MGjRI6GwRQroIlRQb+PLLLzFv3jy89dZbOHnyJPr164eUlJRWp30RQhwHlRQbSEhIwH333ccPPjWZTAgPD8eLL76I1157TeDcEUK6Ag3erqPX65GdnY20tDR+m1QqRVJSksUMD45Op7NYvNVkMqGoqAh+fn42WxqeEGI9xhjKy8sRGhra4px6Cop1bt26BaPRiKCgIIvtQUFB/NzNhpYvX2714geEEOH873//41c/bwoFRSulpaVZrF1YWlqKbt264cMPP0Rqaqr4nl9RXQrFhnsAAIbZuYC6+el4YmYwGHDo0CGMGDFCfPfYDplMDA+vPYqS6lp8PKk/7u3m7bD3uLy8HJGRkS1ORQUoKPL8/f0hk8kaPTK0sLAQwcHBjdKrVKomn2Lm6uoKPz8/8f0xXf4VUEkAr3AgrLvQubGawWAQ7z22U0P7RGBX7nXk3jIi6V4/h73H3LW01rxFvc91lEolBgwYYDG9zGQy4cCBA0hMTBQwZ53kxhnza2DvltMRpzPkbvNSYcfOi391n65AJcUG5s2bh8mTJ2PgwIEYNGgQ3n33XVRWVmLq1KlCZ63jCs2PEkBgrLD5IKIz+G7zCuSnr5airMYAF+GemioKFBQb+MMf/oCbN2/izTffhFarxT333IO9e/c26nyxSzd+Mb8GxQmbDyI6Yd4u6BHghks3K5F18TYe6mnd4wgcBVWf7zB79mxcvnwZOp0OP/zwQ5PPLLE7jFH1mbRoSDRXhb4lcE6ER0HRGZRfB2pKAYkM8G/5+SPEOXFV6KPUrkhB0Slw7Yl+0YC8cY85Iff38IVMKkHB7SpcLa4WOjuCoqDoDLiqcxB1spCmeagVuDfcGwDwn4u3W0nt2CgoOgOukyWQOllI8wbXDc35zwUKisTRFVInC2ndkLp2xaxLRTA58TIxFBQdnbEWuFn3LGeqPpMW9LvLCx5qOUqqDbhaKXRuhENB0dEV5wNGHaBwBby7C50bImJymRSJPcxjFM+WOO9KTxQUHR1XdQ7oBbSwXBIhQP2Uv3OlFBSJo+JnslDVmbSOG6+YXy5Blb5W4NwIg4Kio+NnslBQJK3r7ueKMG81jEyCEwXFQmdHEBQUHR0tBEHaQSKRYHC0uV3xPxeLBM6NMCgoOjJDNVB0yfyegiJpowej6oKik45XpKDoyG6eBcAAVz/APVDo3BA7cX8PX0jA8OuNChSW1QidnS5HQdGRNaw608O0SBv5uCpxl5v5vTOumkNB0ZHdoPZEYp0Yb/OUlmMXKCgSR8IFRRqOQ9qpl1d9UHS2R8NTUHRkty6YXwN6CZsPYnciPRhUcilulutwpahK6Ox0KQqKjooxoEJrfu8ZKmxeiN2RS4FIP1cAwKWbzjURmoKio6ouBox683t3B3jGDOlyPQLMvS0Xb1YInJOuRUHRUZVfN7+6+NJq28Qqkf4UFIkjKa+rOnuECJsPYrd68EGRqs/EEfBBkarOxDpcULxEJUXiELjqM5UUiZUi/c0dLbcq9CitMgicm65DQdFRVRSaX6mThVjJTSVHsKcaAHDxlvOUFikoOioqKZJOEBXIVaGdp12RgqKjKq8rKXoEC5sPYtd6+LsDcK4eaAqKjorvaKGgSKwXFeB8nS0UFB1Rw9ksFBRJB/QI4EqKVH0m9oxms5BOEhVoDoqXb1ei1mgSODddg4KiI+KqzjSbhXRQiKcaaoUUBiPD/4qrhc5Ol6Cg6Ij4nmeqOpOOkUol9Z0tN5yjXZGCoiOqoJ5n0nm4hSEuOclYRQqKjogrKbpTUCQdF8V1ttxwjs4WCoqOiMYokk5EJUVi/6hNkXSiKCcblkNB0RFRmyLpRFxJsahSj+JKvcC5sT0Kio6I2hRJJ3JVyhHqZV4Ywhmq0BQUHQ1j1KZIOp0zzWyhoOhoqosBo878nmazkE4S5UTPa6Gg6Gj4Z7P4AAq1sHkhDoOb7ne+kIIisTcl/zO/eoULmw/iUHqHeAIAfv6tTOCc2B4FRUdTcsX86t1N2HwQh9I7xBMSCaAtq8GtCp3Q2bEpCoqOpuSy+dU7Qth8EIfirpIj0s/crnjGwUuLFBQdTWld9ZlKiqSTxYaaq9BnfisVOCe2RUHR0VD1mdhInzAvAMCZa1RSJPaEgiKxkTgqKRK7o6sAqm6b33tT7zPpXHGh5pJiwe0qlNc47nOgRR0Uly9fjvvuuw8eHh4IDAzEuHHjcO7cOYs0w4cPh0Qisfh54YUXLNJcuXIFqampcHV1RWBgIBYsWIDa2lqLNIcPH0b//v2hUqkQHR2N9PR0W19e5+PaE9XegNpL2LwQh+PrpuSn+zny0BxRB8UjR45g1qxZOH78ODQaDQwGA5KTk1FZaTnVaMaMGbh+/Tr/s3LlSn6f0WhEamoq9Ho9MjMz8cknnyA9PR1vvvkmnyY/Px+pqakYMWIEcnJyMGfOHEyfPh379u3rsmvtFHzVmUqJxDZi60qLjtwDLRc6Ay3Zu3evxef09HQEBgYiOzsbQ4cO5be7uroiOLjpeb7ff/89fv75Z+zfvx9BQUG45557sGzZMixcuBCLFy+GUqnEpk2bEBkZidWrVwMAevfujWPHjmHt2rVISUmx3QV2Nj4o0nAcYht9wjyx/5dC5Dlwu6Kog+KdSkvNvwhfX1+L7Vu3bsVnn32G4OBgjB07FosWLYKrqysAICsrC/Hx8QgKqp8HnJKSgj/+8Y84c+YM7r33XmRlZSEpKcnimCkpKZgzZ06zedHpdNDp6gexlpXV/89pMAjT3iItyocMgNEzDCaB8mBr3L0V6h47g5buca/AurGK10rt7nfQ1vzaTVA0mUyYM2cOHnzwQfTp04ffPnHiRERERCA0NBS5ublYuHAhzp07h6+//hoAoNVqLQIiAP6zVqttMU1ZWRmqq6vh4uLSKD/Lly/HkiVLmsyrRqOx/kI7YGD+CYQB+PlaBS7t2SNIHrqKUPfYmTR1j0t0ACDH+cJyfLtrDxSiboCzVFVV1aZ0dhMUZ82ahby8PBw7dsxi+8yZM/n38fHxCAkJwcMPP4yLFy8iKirKZvlJS0vDvHnz+M9lZWUIDze35Y0cORIKhcJm526O7OO1QAnQOzEFvWLGdPn5u4LBYIBGoxHsHjuDlu4xYwzvnj2M4ioDIu95EH3vsp8OvYa1uZbYRVCcPXs2du3ahYyMDNx1110tpk1ISAAAXLhwAVFRUQgODsaPP/5okaaw0LzeINcOGRwczG9rmMbT07PJUiIAqFQqqFRNP1NZoVAI8w+2rvdZ7t8DcPCAIdg9diLN3eM+YV44ev4Wzt2owoBIfwFyZp22/r2IuvDLGMPs2bPxzTff4ODBg4iMjGz1Ozk5OQCAkJAQAEBiYiL++9//4saNG3wajUYDT09PxMbG8mkOHDhgcRyNRoPExMTOuhTb01cCVbfM72mFHGJD3HhFR+1sEXVQnDVrFj777DNs27YNHh4e0Gq10Gq1qK6uBgBcvHgRy5YtQ3Z2NgoKCvDdd99h0qRJGDp0KPr27QsASE5ORmxsLJ599lmcPn0a+/btwxtvvIFZs2bxJb0XXngBly5dwquvvoqzZ89i48aN2L59O+bOnSvYtbcbt2SY2gtw8RY2L8Sh9QnjZrY45rAcUQfF999/H6WlpRg+fDhCQkL4ny+//BIAoFQqsX//fiQnJ6NXr16YP38+Hn/8cfzrX//ijyGTybBr1y7IZDIkJibimWeewaRJk7B06VI+TWRkJHbv3g2NRoN+/fph9erV+PDDD+10OA5N7yO2xZUUz14vQ63RJHBuOp+o2xQZYy3uDw8Px5EjR1o9TkREBPa00hs7fPhwnDp1ql35ExVuyTAvCorEtiJ8XeGukqNCV4uLNysRE+whdJY6lahLiqQdqKRIuohUKkFsiOMuDkFB0VFQUCRdiFtbMc8BlxGjoOgoaHFZ0oX4tRWppEhEi0qKpAtxayv+/FsZTKaW2/7tDQVFR6CvAipvmt/TCjmkC0QHukMpl6JcV4v/Fbdt+py9oKDoCLiqs8rL/LxnQmxMIZOiV12vs6ONV6Sg6Aio6kwEEMd3tjhWuyIFRUfAP9aUqs6k68Q56IKzFBQdAZUUiQAaPsiqtYkW9oSCoiOgoEgE0DvEEzKpBLcq9LhRrmv9C3aCgqIjKKExiqTrqRUyRAXUrcTtQOMVKSg6AiopEoHwy4g50MwWCor2zlANVNatFUlBkXSxhu2KjoKCor3jqs4qT/PzngnpQlRSJOLDVZ29wgGJRNi8EKfDLQxxraQaJVV6gXPTOSgo2jt+jCJVnUnX83JRoJuv+XHCjjJekYKivaNOFiIwR2tXpKBo72jJMCKw+mXEqKRIxIBKikRgsQ42B5qCor2joEgE1qeuB/rSrUpU6WsFzk3HUVC0Z4ZqoKLQ/J6CIhFIgIcKgR4qMAb8ct3+q9AUFO1Z6VXzq9Kd1lEkgqrvbKGgSITUcDgOjVEkAuI6WxyhXZGCoj2j9kQiElRSJOJAQZGIBDfd79fCcuhrTQLnpmMoKNozWjKMiMRdPi7wVMthMDL8WlgudHY6hIKiPaOSIhEJiUTClxZ/tvMqNAVFe8bNZvGiZ7MQ4fHPgrbzYTkUFO0VY/XPenYPFDYvhADo7m9ehfuqnT8HmoKivdKVA6a62QMuvsLmhRAAYd4uAIDfSmoEzknHUFC0V9VF5le5GlC6CpsXQgCEeKsBAL+VVguck46hoGivquqCIpUSiUiE1pUUS6oMdj0HmoKiveJKiq4UFIk4eKoV8FDJAdh3FZqCor2qKja/0pxnIiJ8FbrEfqvQFBTtFV9S9BM2H4Q0wFWhr9txuyIFRXtVddv8StVnIiJcULxG1WfS5aijhYhQqBdVn4lQqKOFiBBVn4lwqKRIRCjUAQZwU1C0V1RSJCIU6sUFxWowxgTOjXUoKNorfkgOBUUiHkFeKkgkgK7WhKJKvdDZsQoFRXtFJUUiQiq5DP7uKgD2W4WmoGiPanWAvsL8noIiERm+XdFOO1soKNojrpNFIgVUXsLmhZA7hNn5rBYKivaIqzq7+ABS+hUScQlp0Nlij+hflD2i4ThExOqrz9Sm2OkWL14MiURi8dOrVy9+f01NDWbNmgU/Pz+4u7vj8ccfR2FhocUxrly5gtTUVLi6uiIwMBALFixAba3lskaHDx9G//79oVKpEB0djfT09K64POtRJwsRMao+21hcXByuX7/O/xw7dozfN3fuXPzrX//Cjh07cOTIEfz222947LHH+P1GoxGpqanQ6/XIzMzEJ598gvT0dLz55pt8mvz8fKSmpmLEiBHIycnBnDlzMH36dOzbt69Lr7NdqKRIRIyrPl+3095nudAZaI1cLkdwcHCj7aWlpfjoo4+wbds2PPTQQwCALVu2oHfv3jh+/Djuv/9+fP/99/j555+xf/9+BAUF4Z577sGyZcuwcOFCLF68GEqlEps2bUJkZCRWr14NAOjduzeOHTuGtWvXIiUlpUuvtc2opEhEjKs+F5bXwGA0QSETfdnLguhze/78eYSGhqJHjx54+umnceWK+bGe2dnZMBgMSEpK4tP26tUL3bp1Q1ZWFgAgKysL8fHxCAoK4tOkpKSgrKwMZ86c4dM0PAaXhjuGKFVRUCTi5eemhFIuBWOA1g7bFUVdUkxISEB6ejpiYmJw/fp1LFmyBEOGDEFeXh60Wi2USiW8vb0tvhMUFAStVgsA0Gq1FgGR28/taylNWVkZqqur4eLi0mTedDoddDod/7msrP6xjgaDwcorbhtZ5S1IARhV3jDZ+Fxiw91bW99jZ9YZ9zjYU4UrRdW4fKscwdio6r0AACAASURBVB6Kzspah7T1ekQdFEePHs2/79u3LxISEhAREYHt27c3G6y6yvLly7FkyZIm92k0Gpue+4FLuQgAkHPpBq6W7LHpucTK1veYdOweezApACm+OfgDboeIYw50VVXbHr0q6qB4J29vb/Ts2RMXLlzAyJEjodfrUVJSYlFaLCws5Nsgg4OD8eOPP1ocg+udbpjmzh7rwsJCeHp6thh409LSMG/ePP5zWVkZwsPND6UfOXIkFArb/e8o32DuKOo39BH07ZZos/OIkcFggEajsfk9dmadcY9/VV3AmcOXAN9wjBnTp5NzaJ2GtbmW2FVQrKiowMWLF/Hss89iwIABUCgUOHDgAB5//HEAwLlz53DlyhUkJpoDRWJiIt5++23cuHEDgYHmB8ZrNBp4enoiNjaWT7Nnj2VpS6PR8MdojkqlgkqlanKfQqGw3T9YkxEouwoAkPv3AJw0MNj0HhMAHbvHfcPNzw765XqFaH5Pbc2HqDtaXnnlFRw5cgQFBQXIzMzEo48+CplMhgkTJsDLywvTpk3DvHnzcOjQIWRnZ2Pq1KlITEzE/fffDwBITk5GbGwsnn32WZw+fRr79u3DG2+8gVmzZvEB7YUXXsClS5fw6quv4uzZs9i4cSO2b9+OuXPnCnnpzSvXAqZaQCoHPEKEzg0hTYoL9QQA/FpYDl2tUeDctI+oS4pXr17FhAkTcPv2bQQEBGDw4ME4fvw4AgICAABr166FVCrF448/Dp1Oh5SUFGzcuJH/vkwmw65du/DHP/4RiYmJcHNzw+TJk7F06VI+TWRkJHbv3o25c+di3bp1uOuuu/Dhhx+KdzhOibn3HV53AVKZsHkhpBlh3i7wdlWgpMqA84UV6BNmP3P0RR0Uv/jiixb3q9VqbNiwARs2bGg2TURERKPq8Z2GDx+OU6dOWZXHLscHxXBh80FICyQSCeJCPfGfC7eRd63UroKiqKvPpAlcUPSOEDYfhLQiLtQcCM/81rYODrGgoGhvSi6bX727CZsPQlrBtSue+a1U4Jy0DwVFe8OXFCkoEnHjSoq/XC+H0SSOsYptQUHR3lBQJHYi0t8NrkoZqg1G5N+qEDo7bUZB0Z6YTECpeYwiBUUidjKpBL1DzFXovGv2065IQdGeVGgBk4HGKBK7YY/tihQU7QlXdfYMBWSiHk1FCACgT127IpUUiW3QcBxiZ2IblBQZs4/OFgqK9oSG4xA70zPIAwqZBGU1tbhabB+PJ6CgaE+o55nYGaVcip5BHgDsp12RgqI9oaBI7FB9Z4t9tCtSULQnFBSJHeLmPeddo5Ii6Uw0RpHYKSopEtuoKASMekAiAzxChc4NIW3WO8QTEglwo1yHG+Xif5AVBUV7wY9RDKMxisSuuCrl6OHvBsA+SosUFO0FtScSO8YvI2YH7YoUFO0FP0aRFpcl9qdPmP20K1JQtBeVN82vHsHC5oMQK0QFuAMArhS17TGjQqKgaC+qisyvLr7C5oMQK4R6mx8XfL2UOlpIZ6muC4quFBSJ/Qn1MgfFoko9qvXifrofBUV7QSVFYsc8XeRwU5qfPvlbqbjnQFNQtBdUUiR2TCKR1FehS8RdhaagaC+4kqKrn7D5IMRKIXVB8bcSKimSjjIaAF3dUAaqPhM7FeatBgBco6BIOqy6uO6NBHDxFjQrhFiL62y5Tm2KpMO4qrPaC5DKhM0LIVaqrz5TmyLpKOpkIQ4gtK76TL3PpONoOA5xAGENOlrE/LwWCor2gEqKxAEEe5lLijUGE4qrDALnpnkUFO1B1W3zK5UUiR1TyWXwd1cBEPewHAqK9oDGKBIHwbcrUlAkHcJXn32EzQchHcQNy6GgSDqmqm6cIlWfiZ2zh9VyKCjaA+poIQ4i1A5mtVBQtAc0JIc4CCopks5BJUXiIELtYFEICopixxiVFInDCK0bq1hYVoNao0ng3DSNgqLY1ZQCrG6lYiopEjvn766CQiaBiQGF5Tqhs9MkCopix1WdFa6AwkXYvBDSQVKpBCEiH5ZDQVHsaDgOcTAhXuIewE1BUex0dQ8PV3sJmw9COkmQpzko3qTqM7GKrsL8qnIXNh+EdBIfVwUAoLhKL3BOmkZBUez0dUFR6SZsPgjpJD5uSgAQ7Uo5FBTFjispKqmkSByDLxcUK6mkSKzBlRRVHsLmg5BO4uNqDopFFBSJVfRUUiSOhS8pUpsisQp1tBAH4813tFCbolW6d+8OiUTS6GfWrFkAgOHDhzfa98ILL1gc48qVK0hNTYWrqysCAwOxYMEC1NbWWqQ5fPgw+vfvD5VKhejoaKSnp3fVJbZMX2l+pY4W4iAatimK8VktcqEz0JoTJ07AaDTyn/Py8jBy5EiMHz+e3zZjxgwsXbqU/+zq6sq/NxqNSE1NRXBwMDIzM3H9+nVMmjQJCoUCf/3rXwEA+fn5SE1NxQsvvICtW7fiwIEDmD59OkJCQpCSktIFV9kCfbn5VUltisQxcG2KtSaGcl0tPNUKgXNkSfRBMSAgwOLzihUrEBUVhWHDhvHbXF1dERwc3OT3v//+e/z888/Yv38/goKCcM8992DZsmVYuHAhFi9eDKVSiU2bNiEyMhKrV68GAPTu3RvHjh3D2rVrhQ+KVH0mDkatkMFVKUOV3ojiSr3ogqLoq88N6fV6fPbZZ3juuecgkUj47Vu3boW/vz/69OmDtLQ0VFVV8fuysrIQHx+PoKAgfltKSgrKyspw5swZPk1SUpLFuVJSUpCVlWXjK2oD6mghDogrLYqxXVH0JcWGdu7ciZKSEkyZMoXfNnHiRERERCA0NBS5ublYuHAhzp07h6+//hoAoNVqLQIiAP6zVqttMU1ZWRmqq6vh4tJ4IQadTgedrn6aUllZGf/eYOi8X7RcVwEJgFqZGqwTj2uvuHvbmfeYWOqKe+ztKse1EuBmWRUMhq5pL2/r9dhVUPzoo48wevRohIaG8ttmzpzJv4+Pj0dISAgefvhhXLx4EVFRUTbLy/Lly7FkyZIm92k0mk47T1LJDbgByMzORfFZcU6gF0Jn3mPSNFveY2OlFIAUh7N+QtWFrulsaViDbIndBMXLly9j//79fAmwOQkJCQCACxcuICoqCsHBwfjxxx8t0hQWFgIA3w4ZHBzMb2uYxtPTs8lSIgCkpaVh3rx5/OeysjKEh4cDAEaOHAmFonPaSeTn5gJ6IHFYMhDYu1OOac8MBgM0Gk2n3mNiqSvu8f7KXJzN1aJbdG+MebC7Tc5xp4a1uZbYTVDcsmULAgMDkZqa2mK6nJwcAEBISAgAIDExEW+//TZu3LiBwMBAAOb/AT09PREbG8un2bNnj8VxNBoNEhMTmz2PSqWCSqVqcp9Coei8P6a6jhaFmzdAQYDXqfeYNMmW99jP3bxSTmmNsct+j209j110tJhMJmzZsgWTJ0+GXF4fxy9evIhly5YhOzsbBQUF+O677zBp0iQMHToUffv2BQAkJycjNjYWzz77LE6fPo19+/bhjTfewKxZs/ig9sILL+DSpUt49dVXcfbsWWzcuBHbt2/H3LlzBblentEAGOvaLamjhTgQMXe02EVQ3L9/P65cuYLnnnvOYrtSqcT+/fuRnJyMXr16Yf78+Xj88cfxr3/9i08jk8mwa9cuyGQyJCYm4plnnsGkSZMsxjVGRkZi9+7d0Gg06NevH1avXo0PP/xQ+OE4XM8zQEGROBRft7pZLSKc/2wX1efk5OQmR76Hh4fjyJEjrX4/IiKiUfX4TsOHD8epU6eszqNNcGMUZUpArhQ2L4R0Im75sCIRzn+2i5Ki06IxisRB+bqKd/kwCopixs97pqBIHIs3tSkSq+jq5j3TFD/iYBouHya2RSEoKIoZVZ+Jg+KWDzOaGMpqaltJ3bUoKIoZV32mkiJxMGqFDG5KGQDxtStSUBSzqiLzq9pb2HwQYgPco07F9vxnCopiVn7d/OoRImw+CLGBHgHmhSAu3qxoJWXXoqAoZhV187E9glpOR4gd6hFgbha6eLNS4JxYoqAoZlRSJA4sikqKpN3K60qK7lRSJI6HKyleopIiabNy8yK4VFIkjiiqLiheK6lGtd7YSuquQ0FRrPRVgK7U/J7aFIkD8nVT8uMVL90STxWagqJYVdSVEhWugMpT2LwQYiNRIqxCU1AUq4btiQ0e0kWII+nhL77OFgqKYkU9z8QJRAVSSZG0FY1RJE6ASoqk7aikSJxAw5KiySSO1XIoKIoVjVEkTqCbryvkUgmqDUZoy2qEzg4ACoriRSVF4gQUMim6+bkCEE8VmoKiWFGbInESYhuWQ0FRrKikSJyE2FbLoaAoRoZqoKZuNgu1KRIHRyVF0jpuzrPcBVB7CZsXQmxMbKvlUFAUo4btiTSbhTi4Hv7mkuL10hpU6oR/XgsFRTGi9kTiRHzclPzT/fJvCV+FpqAoRjRGkTgZMVWhKSiKEZUUiZPhqtBieDQBBUUxojGKxMlEBVJJkbSk8pb51dVf2HwQ0kW6+ZpntVwrFv5xpxQUxUhf97+lykPYfBDSRXzdVACAkiq9wDmhoChOfFB0FzYfhHQRXzfzYwmKKikokqbo6oKikkqKxDn4uJqH5JTV1KLWaBI0LxQUxYgrKSrdhM0HIV3Ey0XBvy+pNgiYEwqK4qSj6jNxLnKZlA+MxQJXoeXWftFkMuHChQu4ceMGTCbL4u7QoUM7nDGnZawFaut64Kj6TJyIr5sSpdUGwdsVrQqKx48fx8SJE3H58mUwZrmEuEQigdEongdb2x1Dg8GrVFIkTsTHVYF8AMUC90BbFRRfeOEFDBw4ELt370ZISAgktGhB5+GqzlIFIFcJmxdCuhA3/7m4Stg2RauC4vnz5/HVV18hOjq6s/NDqJOFOCnvuh5ooavPVnW0JCQk4MKFC52dFwI06GSh9kTiXPiSor20Kebm5vLvX3zxRcyfPx9arRbx8fFQKBQWafv27dt5OXQ2fEmR2hOJc+HGKhbZS5viPffcA4lEYtGx8txzz/HvuX3U0dJBNJuFOCluVkuJvbQp5ufn2zIfhKOjNkXinMTSptjmoBgREcG/z8jIwAMPPAC53PLrtbW1yMzMtEhL2klfbn6l6jNxMvW9z3bY0TJixAgUFRU12l5aWooRI0Z0OFNOTV83TpE6WoiT8RFJSdGqoMi1Hd7p9u3bcHOjal+H6KijhTgnrqRYXlMLg4CLQrRrnOJjjz0GwNypMmXKFKhU9YOLjUYjcnNz8cADD3RuDp0NdbQQJ+XlooBEAjBm7mwJ8BBm8kK7gqKXl/kZxIwxeHh4wMXFhd+nVCpx//33Y8aMGZ2bQ2ej49oUqcRNnItMKoGXiwIlVQYUV+kFC4rtqj5v2bIFW7ZswVtvvYWPPvqI/7xlyxZs3rwZaWlp8Pdv+xL6GRkZGDt2LEJDQyGRSLBz506L/YwxvPnmmwgJCYGLiwuSkpJw/vx5izRFRUV4+umn4enpCW9vb0ybNg0VFZbPecjNzcWQIUOgVqsRHh6OlStXNsrLjh070KtXL6jVasTHx2PPnj3tuDOdiGtTpMUgiBPyFUG7olVtim+99Rbc3Nxw48YNHD16FEePHsWNGzfafZzKykr069cPGzZsaHL/ypUr8fe//x2bNm3CDz/8ADc3N6SkpKCmpoZP8/TTT+PMmTPQaDTYtWsXMjIyMHPmTH5/WVkZkpOTERERgezsbKxatQqLFy/GBx98wKfJzMzEhAkTMG3aNJw6dQrjxo3DuHHjkJeX1+5r6jCqPhMn5lPXrijoYwmYFcrKytgzzzzD5HI5k0gkTCKRMLlczp5++mlWUlJizSEZAPbNN9/wn00mEwsODmarVq3it5WUlDCVSsU+//xzxhhjP//8MwPATpw4waf597//zSQSCbt27RpjjLGNGzcyHx8fptPp+DQLFy5kMTEx/Ocnn3ySpaamWuQnISGBPf/8823Of2lpKQPAtm3bxvR6fZu/18hHoxh7y5OxvG9aT+uE9Ho927lzZ8fuMWmRkPd4WvoJFrFwF9t6/HKnH5v7N1paWtpiOqsWhJg+fTpOnTqFXbt2ITExEQCQlZWFl19+Gc8//zy++OKLDgfr/Px8aLVaJCUl8du8vLyQkJCArKwsPPXUU8jKyoK3tzcGDhzIp0lKSoJUKsUPP/yARx99FFlZWRg6dCiUSiWfJiUlBe+88w6Ki4vh4+ODrKwszJs3z+L8KSkpjarzDel0Ouh0Ov5zWVkZ/95gsH5EvlxXDgmAWpkLWAeO46i4e9uRe0xaJuQ99nKRAQBulVd3+vnbejyrguKuXbuwb98+DB48mN+WkpKC//u//8OoUaOsOWQjWq0WABAUZPns46CgIH6fVqtFYGCgxX65XA5fX1+LNJGRkY2Owe3z8fGBVqtt8TxNWb58OZYsWdLkPo1G09rlNevhkhtwB5CVnYuiczWtpndWHbnHpG2EuMfFWikAKU6e+RV7Ks926rGrqqralM6qoOjn58f3RDfk5eUFHx8faw5pd9LS0ixKl2VlZQgPDwcAjBw5stEiGW0l//UVQAfcP2wkEBTXKXl1JAaDARqNpkP3mLRMyHt89Wg+Dv52Hj5BYRgzJr5Tj92wNtcSq4LiG2+8gXnz5uHTTz9FcHAwAHOpa8GCBVi0aJE1h2yEO25hYSFCQkL47YWFhbjnnnv4NHd28NTW1qKoqIj/fnBwMAoLCy3ScJ9bS8Ptb4pKpbIYp9mQQqGw/o+prqNF4eYN0D/6ZnXoHpM2EeIeB3iYh/mVVNd2+rnbejyrep/ff/99HD9+HN26dUN0dDSio6PRrVs3ZGZmYvPmzejfvz//Y63IyEgEBwfjwIED/LaysjL88MMPfDtmYmIiSkpKkJ2dzac5ePAgTCYTEhIS+DQZGRkW7QkajQYxMTF8qTYxMdHiPFwa7jxdxmQEDHVFfBqSQ5yQt2vd858FXCnHqpLiuHHjOuXkFRUVFovV5ufnIycnB76+vujWrRvmzJmDv/zlL7j77rsRGRmJRYsWITQ0lD9/7969MWrUKMyYMQObNm2CwWDA7Nmz8dRTTyE0NBQAMHHiRCxZsgTTpk3DwoULkZeXh3Xr1mHt2rX8eV9++WUMGzYMq1evRmpqKr744gv89NNPFsN2uoS+wfNZaPA2cUKiWGi20/u92+HQoUMMQKOfyZMnM8bMw3IWLVrEgoKCmEqlYg8//DA7d+6cxTFu377NJkyYwNzd3ZmnpyebOnUqKy8vt0hz+vRpNnjwYKZSqVhYWBhbsWJFo7xs376d9ezZkymVShYXF8d2797drmvplCE5pdfMw3GW+DJmMll3DAdHQ3JsT8h7fOFGOYtYuIv1eWtvpx/bpkNyAKCkpARfffUVLl68iAULFsDX1xcnT55EUFAQwsLC2nSM4cOHN3oaYEMSiQRLly7F0qVLm03j6+uLbdu2tXievn374ujRoy2mGT9+PMaPH99yhm2t4WIQ9DAw4oS4GS3cohAKWdc/mt6qoJibm4ukpCR4eXmhoKAAM2bMgK+vL77++mtcuXIF//jHPzo7n86B1lIkTs6zwaIQxVV6BHqouzwPVoXhefPmYcqUKTh//jzU6vpMjxkzBhkZGZ2WOafDr6VIQZE4J5lUAm8Xc2dLcaUwnS1WBcUTJ07g+eefb7Q9LCysxQHPpBW0liIh/PxnoVbgtiooqlSqJgdC/vrrrwgICOhwppwWLQZBCN+uKFQPtFVB8Xe/+x2WLl3Kj/2TSCS4cuUKFi5ciMcff7xTM+hUdNSmSAhXUhTqUadWBcXVq1ejoqICAQEBqK6uxrBhwxAdHQ0PDw+8/fbbnZ1H50HPfCYEPq5cm6IwQdGq3mcvLy9oNBr85z//wenTp1FRUYH+/ftbrGhDrEAdLYTUlxQF6mhpd1A0mUxIT0/H119/jYKCAkgkEn5KHmvmgVakjaijhRC+TVGohWbbVX1mjOF3v/sdpk+fjmvXriE+Ph5xcXG4fPkypkyZgkcffdRW+XQONE6REMHbFNtVUkxPT0dGRgYOHDjQ6PnOBw8exLhx4/CPf/wDkyZN6tRMOg0d9T4T4mNPvc+ff/45Xn/99SYfeP/QQw/htddew9atWzstc06Hf2gVBUXivHzduJVy7CAo5ubmtriy9ujRo3H69OkOZ8pp0ThFQviSYok9zGgpKipqtGx/Q0FBQSguLu5wppwWjVMkhF8+rFxXC32tqcvP366gaDQaIZc33wwpk8lQW1vb4Uw5LRqnSAg81QpI6waxCNED3a6OFsYYpkyZ0uwy/A2fbkesQOMUCYFUKoG3qxJFlXoUVekR6Nm1K+W0KyhOnjy51TTU89wB/DhFWnWbODcfVwWKKvWCrJTTrqC4ZcsWW+WDMAYYqPeZEADwdlUCqERpdddXn7t+WVvSNGOD/xFlSuHyQYgIuCplAIBqg7HLz01BUSyMDf5HlDfdZkuIs1Ar6oKiXuS9z8SGGgZFKikSJ+eioJIi4YKiRApIZcLmhRCBcUGxhoKiE6utG84ko6ozIS5cm6KegqLz4jpaqOpMSH2bIpUUnRhXfZYphM0HISKgVphDE1WfnZmxrvpMPc+EUEcLQYPqM5UUCeHaFKmk6Mz46jO1KRJSP06RgqLzot5nQnhUfSZUfSakgfqgSDNanBdXfaaOFkLq2xSp+uzEaEgOITwap0ioo4WQBmiaH6GOFkIacKGlwwh1tBBSj0qKhKrPhDTABUWDkcFg7NoeaAqKYsFP86OgSIhKUR+aurq0SEFRLGiVHEJ4KrkUkrrHnHZ1uyIFRbGg6jMhPIlEUt+u2MWPJKCgKBZ87zMFRUIA4ab6UVAUCyopEmKBG5ZTqa/t0vNSUBQLfd0zn1X0zGdCAMDb1Tw8rbTK0ErKzkVBUSx05eZXJQVFQgDAx9Vcayqq1LeSsnNRUBQLfYX5lYIiIQAAXzdzUCyuoqDonHR1QZGqz4QAqC8pUlB0VlRSJMRCffWZ2hSdExcUVR7C5oMQkfB1M3e0FFObopPiqs9KN2HzQYhI+NS1KRY5U/U5IyMDY8eORWhoKCQSCXbu3MnvMxgMWLhwIeLj4+Hm5obQ0FBMmjQJv/32m8UxunfvDolEYvGzYsUKizS5ubkYMmQI1Go1wsPDsXLlykZ52bFjB3r16gW1Wo34+Hjs2bPHNhfdHKo+E2LBl2tTdKaSYmVlJfr164cNGzY02ldVVYWTJ09i0aJFOHnyJL7++mucO3cOv/vd7xqlXbp0Ka5fv87/vPjii/y+srIyJCcnIyIiAtnZ2Vi1ahUWL16MDz74gE+TmZmJCRMmYNq0aTh16hTGjRuHcePGIS8vzzYXfqdaff3gbepoIQQA4M13tHRtm6K8S892h9GjR2P06NFN7vPy8oJGo7HYtn79egwaNAhXrlxBt27d+O0eHh4IDg5u8jhbt26FXq/Hxx9/DKVSibi4OOTk5GDNmjWYOXMmAGDdunUYNWoUFixYAABYtmwZNBoN1q9fj02bNnXGpbaMKyUCVFIkpE7DITmMMUi4FSJsTNCg2F6lpaWQSCTw9va22L5ixQosW7YM3bp1w8SJEzF37lzI5eZLy8rKwtChQ6FU1k+fS0lJwTvvvIPi4mL4+PggKysL8+bNszhmSkqKRXX+TjqdDjqdjv9cVlbGvzcY2vk/W1UJFACYTIVaEwBT1/7PaE+4e9vue0zaTCz32L1uvWWjiaGovBqeLh1bgLmt12M3QbGmpgYLFy7EhAkT4OnpyW9/6aWX0L9/f/j6+iIzMxNpaWm4fv061qxZAwDQarWIjIy0OFZQUBC/z8fHB1qtlt/WMI1Wq202P8uXL8eSJUua3HdnCbc1HtVX8RAAvUSJvV3dlmmn2nuPSfuJ4R6rpDLoTBJ8s0eDAJeOHauqqqpN6ewiKBoMBjz55JNgjOH999+32NewhNe3b18olUo8//zzWL58OVQq2z3vJC0tzeLcZWVlCA8PBwCMHDkSCkXb/1eTXD0BnAWU7r4YM2ZMp+fVkRgMBmg0mnbfY9J2YrrHK3/JwLWSGvQb9ADuCfdu/QstaFiba4nogyIXEC9fvoyDBw9alBKbkpCQgNraWhQUFCAmJgbBwcEoLCy0SMN95tohm0vTXDslAKhUqmaDrkKhaN8fk6kGACBReQj+R2gv2n2PSbuJ4R77uqlwraQG5XpTh/PS1u+LepwiFxDPnz+P/fv3w8/Pr9Xv5OTkQCqVIjAwEACQmJiIjIwMi/YEjUaDmJgY+Pj48GkOHDhgcRyNRoPExMROvJoW0BhFQprEj1XswlktgpYUKyoqcOHCBf5zfn4+cnJy4Ovri5CQEDzxxBM4efIkdu3aBaPRyLfx+fr6QqlUIisrCz/88ANGjBgBDw8PZGVlYe7cuXjmmWf4gDdx4kQsWbIE06ZNw8KFC5GXl4d169Zh7dq1/HlffvllDBs2DKtXr0Zqaiq++OIL/PTTTxbDdmyKxigS0iQfVwFmtTABHTp0iAFo9DN58mSWn5/f5D4A7NChQ4wxxrKzs1lCQgLz8vJiarWa9e7dm/31r39lNTU1Fuc5ffo0Gzx4MFOpVCwsLIytWLGiUV62b9/OevbsyZRKJYuLi2O7d+9u17WUlpYyAGzbtm1Mr9e370Yc38zYW56Mffls+77nhPR6Pdu5c2f77zFpMzHd47e+zWMRC3exFf/+pcPH4v6NlpaWtphO0JLi8OHDwRhrdn9L+wCgf//+OH78eKvn6du3L44ePdpimvHjx2P8+PGtHssm9NxaijTvmZCGuLGKJV041U/UbYpOg1bdJqRJ9W2KFBSdC3W0ENKk+jbFrutooaAoBtTRQkiTuEUhunKlHAqKYsA9n4XWUiTEAld97sreZwqKYsC1KVL1mRAL3OrbJdWGVjteOwsFRTHglg2T225aIiH2yFVlfvaz0cSgN5q65JwUFMWAC4oyZcvpCHEyLgoZ/75GT0HRedTWLUEmMyNRCQAAHe1JREFUo5IiIQ0pZFLIpeZ1FKsNxi45JwVFMTDWDTeQ0QIHhNyJKy1SUHQmVH0mpFkqLijqKSg6D2Nd9Zk6WghpxEVpDlNUUnQmVH0mpFlc9VlHQdGJUPWZkGZRm6IzqqWgSEhz1BQUnRCVFAlplouSOlqcC2M0o4WQFnDV5xoqKToJUy3MC4qDOloIaQK1KTobY4PVP6j6TEgj9eMUaZqfc+Cm+AE0zY+QJlBJ0dlwYxQhAaSyFpMS4oy4wdvUpugsGvY8SyTC5oUQEaKOFmdDPc+EtIjGKTobvqRIPc+ENIXGKTobGrhNSIuoo8XZ0BQ/QlqkpjZFJ0MlRUJaRCVFZ0MdLYS0SE2LzDoZ6mghpEXU0eJs9BXmVwU985mQpni5mAsMpdWGVlJ2DgqKQtPVBUWVu7D5IESkfF3N7e2VemOXdLZQUBSavtL8qqSgSEhTPNRyyOoec1pSZfvSIgVFoXHVZyVVnwlpilQqgY+ruQpdXKVvJXUnnM/mZyAt05WbX1UewuaDEBHzqatCF1dSUHR8fEmRqs+ENIcLikVUUnQCXJsidbQQ0iwft7rqM5UUnYCOSoqEtMbXra6kWEkdLY5PX9emSEGRkGbxbYpUfXYCNE6RkFZRUHQm1NFCSKt8+OozBUXHRx0thLTK143GKToP6mghpFX14xSpo8WxMUYdLYS0Adf7TCVFR2eoBljdA76p+kxIs7zrSopVXbAoBAVFIXHtiQAtHUZICzwbLAph69IiBUUhcVVnhRsgpV8FIc2RSCT1U/1s3ANN/xKFRGMUCWkzrgfa1suHCRoUMzIyMHbsWISGhkIikWDnzp0W+6dMmQKJRGLxM2rUKIs0RUVFePrpp+Hp6Qlvb29MmzYNFRUVFmlyc3MxZMgQqNVqhIeHY+XKlY3ysmPHDvTq1QtqtRrx8fHYs2dP51/wnWiMIiFt5hQlxcrKSvTr1w8bNmxoNs2oUaNw/fp1/ufzzz+32P/000/jzJkz0Gg02LVrFzIyMjBz5kx+f1lZGZKTkxEREYHs7GysWrUKixcvxgcffMCnyczMxIQJEzBt2jScOnUK48aNw7hx45CXl9f5F90QlRQJabOumtUit+nRWzF69GiMHj26xTQqlQrBwcFN7vvll1+wd+9enDhxAgMHDgQAvPfeexgzZgz+9re/ITQ0FFu3boVer8fHH38MpVKJuLg45OTkYM2aNXzwXLduHUaNGoUFCxYAAJYtWwaNRoP169dj06ZNnXjFd+BLirSWIiGt6apZLaJvUzx8+DACAwMRExODP/7xj7h9+za/LysrC97e3nxABICkpCRIpVL88MMPfJqhQ4dCqax/rnJKSgrOnTuH4uJiPk1SUpLFeVNSUpCVlWXLS6NVtwlph65qUxS0pNiaUaNG4bHHHkNkZCQuXryI119/HaNHj0ZWVhZkMhm0Wi0CAwMtviOXy+Hr6wutVgsA0Gq1iIyMtEgTFBTE7/Px8YFWq+W3NUzDHaMpOp0OOp2O/1xWVsa/Nxja9kuTVpdBBsCkcIWxjd9xdty9bes9Ju0n1nvsqTI/6vRWeY1VeWvrd0QdFJ966in+fXx8PPr27YuoqCgcPnwYDz/8sIA5A5YvX44lS5Y0uU+j0bTpGHdrcxAL4H/aW8jpio4dB9LWe0ysJ7Z7nF8oASBDwdXfsGfP1XZ/v6qqqk3pRB0U79SjRw/4+/vjwoULePjhhxEcHIwbN25YpKmtrUVRURHfDhkcHIzCwkKLNNzn1tI015YJAGlpaZg3bx7/uaysDOHh4QCAkSNHQqFo/eH20qNngOvAXRE9EDpmTKvpifl/e41G0+Z7TNpPrPe49vR1fHHpv/D08ceYMQNb/8IdGtbmWmJXQfHq1au4ffs2QkJCAACJiYkoKSlBdnY2BgwYAAA4ePAgTCYTEhIS+DR//vOfYTAY+F+wRqNBTEwMfHx8+DQHDhzAnDlz+HNpNBokJiY2mxeVSgWVStXkPoVC0bY/JmaeriRTqCET0R+fPWjzPSZWE9s9dncx9wvU1JqsyldbvyNoR0tFRQVycnKQk5MDAMjPz0dOTg6uXLmCiooKLFiwAMePH0dBQQEOHDiA3//+94iOjkZKSgoAoHfv3hg1ahRmzJiBH3/8Ef/5z38we/ZsPPXUUwgNDQUATJw4EUqlEtOmTcOZM2fw5ZdfYt26dRalvJdffhl79+7F6tWrcfbsWSxevBg//fQTZs+ebdsbYKxrk5QrW05HCIGLwtymWG0w2fZETECHDh1iABr9TJ48mVVVVbHk5GQWEBDAFAoFi4iIYDNmzGBardbiGLdv32YTJkxg7u7uzNPTk02dOpWVl5dbpDl9+jQbPHgwU6lULCwsjK1YsaJRXrZv38569uzJlEoli4uLY7t3727XtZSWljIAbNu2bUyv17ftS3teZewtT8b2L2nXuZyZXq9nO3fubPs9Ju0m1nv8w6XbLGLhLjZ81SGrvs/9Gy0tLW0xnaDV5+HDh4Mx1uz+ffv2tXoMX19fbNu2rcU0ffv2xdGjR1tMM378eIwfP77V83UqY914K1nT1XBCSD2+pKinVXIcVy0XFMXTbkOIWLkozeGqppaCouPiS4rUpkhIa9RUUnQCXFCUU/WZkNZw1WddrQkmU/PNbh1FQVFIRqo+E9JWLkoZ/96WVWgKikKi6jMhbaaW1wdFW1ahKSgKiYIiIW0mlUqgkptDVrUNn9NCQVFItRQUCWkPrgpty4dXUVAUEpUUCWkXrgpdY8NZLRQUhcT3PlNQJKQtuJIiVZ8dFZUUCWmXrhirSEFRSDTNj5B2cVFQR4tjo2l+hLQLdbQ4Oqo+E9IuXbEoBAVFIRnrnhlB0/wIaRO+TZFKig6KW2SWqs+EtIkLBUUHR9VnQtqFKynSOEVHZDICrO4XS0GRkDbhxynqa212DgqKQqmtf2Y0BUVC2oZKio6MqzoDFBQJaSNqU3RkFkGROloIaQsavO3IGnaySCTC5oUQO8EP3qZxig6IpvgR0m40TtGR0RQ/QtrNRUHT/BwXjVEkpN3qS4rU++x4+Cl+FBQJaStaEMKR8VP8KCgS0la0IIQjo+ozIe1GHS2OjKs+U1AkpM3ocQSOrJaqz4S0F1d91teaYDQxm5yDgqJQqPpMSLtxQRGwXWcLBUWh0JP8CGk3lbw+ZFFQdDRUUiSk3aRSCR8YbdWuSEFRKBQUCbGKrccqUlAUSi0FRUKsUT9W0TazWigoCoVKioRYxdZrKlJQFApN8yPEKrYewE1BUSg0zY8Qq9Q/p4WComOh6jMhVrH18mEUFIVC0/wIsYq67pEEFBQdDU3zI8QqXJtiFVWfHYyRVt4mxBp89bmWgqJj4af50TNaCGkPWz+8ioKiUKijhRCr0DhFR0UdLYRYhcYpOirqaCHEKvXjFGman2Oh6jMhVqFxio6KpvkRYhWHblPMyMjA2LFjERoaColEgp07d1rsl0gkTf6sWrWKT9O9e/dG+1esWGFxnNzcXAwZMgRqtRrh4eFYuXJlo7zs2LEDvXr1glqtRnx8PPbs2WObi+bQND9CrKJ25Gl+lZWV6NevHzZs2NDk/uvXr1v8fPzxx5BIJHj88cct0i1dutQi3YsvvsjvKysrQ3JyMiIiIpCdnY1Vq1Zh8eLF+OCDD/g0mZmZmDBhAqZNm4ZTp05h3LhxGDduHPLy8mxz4QCNUyTESrYepyi3yVHbaPTo0Rg9enSz+4ODgy0+f/vttxgxYgR69Ohhsd3Dw6NRWs7WrVuh1+vx8ccfQ6lUIi4uDjk5OVizZg1mzpwJAFi3bh1GjRqFBQsWAACWLVsGjUaD9evXY9OmTR25xObxvc80TpGQ9uCm+TlkSbE9CgsLsXv3bkybNq3RvhUrVsDPzw/33nsvVq1ahf9v795jmjr/P4C/ezulnQVBhcJGEYdzbk6cqNjt91P4BUG/ZAm6MOMlcRfnxooTWbygxssSR7LE4TIvWzan/jEH8w91M0ZkKOBi3YaGZW6MaLyRQcErIJden98ffHts5VZ6oe3p55WQ0PZzTp/z0H445zw3i8XCv6bX6zFnzhxw3OPL1KysLDQ0NODBgwd8TEZGhtM+s7KyoNfrfXQ0oNZnQtzk64YWv54pDsfhw4ehUqmwaNEip+c//PBDTJ8+HVFRUbhw4QKKiorQ3NyMzz77DABgMBiQmJjotE1MTAz/WmRkJAwGA/+cY4zBYBiwPEajEUajkX/c3t7O/242m4c8HqnVBBEACxOBuRBPetnr1pU6Ju4J9DqWinqXNu02WYdVRldjgyYpfvvtt1i2bBnCwsKcni8sLOR/nzp1KjiOw3vvvYfi4mLI5b67NC0uLsaOHTv6fa2iomLI7bO6OhAGoEb/GzoUzV4unfC5UsfEM4Faxy3dACBFe1fPsBpEu7q6XIoLiqR4/vx5NDQ0oKysbMjY1NRUWCwW3Lx5E5MmTYJarUZLS4tTjP2x/T7kQDED3acEgKKiIqeE3N7ejvj4eADAvHnzIJMN3oAirRcBFuB/0/4PGDNxyOMivcxmMyoqKlyqY+KeQK/j5rYefFJXAwvE+M9/slzezvFqbjBBkRQPHDiAlJQUJCcnDxlbV1cHsViM6OhoAIBWq8XmzZthNpv5P3BFRQUmTZqEyMhIPqayshIFBQX8fioqKqDVagd8H7lcPuCZqEwmG/rD9N+GFpn8KSAAP3iBzqU6Jh4J1DpWKXovn81WBpFYAqnEtaYRV4/Fr0nx0aNHuHbtGv/4xo0bqKurQ1RUFDQaDYDe7H706FHs2rWrz/Z6vR6//vor0tPToVKpoNfrsXbtWixfvpxPeEuXLsWOHTvwzjvvYMOGDbhy5Qo+//xzlJSU8PtZs2YN5s6di127diE7OxulpaWora116rbjdSHa0GK1Du8+0JPMZjOkUil6enpgtfrmRnuwkslkkEgk/i6Gz9mH+QFAj8WGUS4mRVf5NSnW1tYiPT2df2y/HF2xYgUOHToEACgtLQVjDEuWLOmzvVwuR2lpKbZv3w6j0YjExESsXbvW6bI2IiICZ86cgU6nQ0pKCsaOHYutW7fy3XEA4JVXXsGRI0ewZcsWbNq0CRMnTsTx48cxZcoU3xy4zQqw/36hQyQpMsZgMBjw8OFDj/ejVqvR2NgIkUjkpdIJx+jRo6FWqwVdN3Lp4yTYbbJilNy7acyvSTEtLQ2MsUFjVq1a5ZTAHE2fPh0XL14c8n2mTp2K8+fPDxqTm5uL3NzcIfflFVaHM6UQGeZnT4jR0dFQKpVuf2ltNhsePXqEUaNGQSwOmh5lPscYQ1dXF1pbWwEAsbGxfi6R74hEIihkEnSbrT7plhMU9xQFx/q4K08onClarVY+IY4ZM8ajfdlsNphMJoSFhVFSfIJCoQAAtLa2Ijo6WtCX0gquNyn6Yvwzfar8wfFMURx4N7K9zX4PUalU+rkkwmev40DtY+gtvuzATUnRH+zjnsUyIITOdoR8nytQhEod+3KoX+h8IwNJiLY8E+Itvpx9m5KiP/CTQQj/0pkQX6DLZ6GhuRSDxptvvsnP0ymTyZCYmIj169ejp6eHj+lvLlD7tjk5OU7PNTY24u2330ZcXBw4jkNCQgLWrFmDe/fuOcWlpaU5zREaExOD3Nxc3Lp1yzcHGmTsfRV9sfYzJUV/oOVNg8r8+fPR3NyM69evo6SkBF999RW2bds27P1cv34dM2bMwNWrV/H999/j2rVr+PLLL1FZWQmtVov79+87xb/77rtobm5GU1MTTpw4gcbGRixfvtxbhxXUwvgzRe+v00JJ0R9oJb+gIpfLoVarER8fj5ycHGRkZLg1WYJOpwPHcThz5gzmzp0LjUaDBQsW4Oeff8a///6LzZs3O8UrlUqo1WrExsZi9uzZyM/Px+XLl711WEHNl0sSUD9Ff6CGFjDG3PpA22w2dJuskJosbvVTVMgkHrXQXrlyBRcuXEBCQsKwtrt//z7Ky8uxc+dOvj+hnVqtxrJly1BWVoZ9+/b1W7779+/jhx9+QGpqqttlFxJf3lOkpOgP1NCCbrMVL2wtH/H3/fvjLCi54X3sT548iVGjRsFiscBoNEIsFmPPnj1OMUuWLOnTWdpoNCI7OxsAcPXqVTDGMHny5H7fY/LkyXjw4AHu3LnDT2ayb98+fPPNN/xoleeeew7l5SNfZ4FI4cN1Wigp+oO9oYXuKQaF9PR07N+/H52dnSgpKYFUKu2zTlBJSUmf2ds3bNjQZ9KKoYa1Olq2bBl/Sd3S0oJPPvkEmZmZuHTpElQqlZtHIwxhdKYoMLTmMxQyCf7+2PW58OxsNhs62jugCle5ffk8XE899RSSkpIA9E52nJycjAMHDjgtjaFWq/kYO5VKxU+AkZSUBJFIhPr6eixcuLDPe9TX1yMyMhLjxo3jn4uIiOD3mZSUhAMHDiA2NhZlZWVYuXLlsI9DSPjO29QlRyCooQUikQhKTurWj4KTuL2tpyM+xGIxNm3ahC1btqC7u9vl7caMGYN58+Zh3759fbYzGAz47rvvsHjx4kHLZ788H877CpUvG1ooKfoDNbQEtdzcXEgkkgGX5h3Inj17YDQakZWVhZqaGjQ2NuL06dOYN28enn76aezcudMpvqurCwaDAQaDAX/88Qfy8vIQFhaGzMxMbx5OULLfU6TO20JBaz4HNalUivz8fHz66afo7Ox0ebuJEyeitrYWEyZMwBtvvIFnn30Wq1atQnp6OvR6PaKiopziv/76a8TGxiI2Nhbp6em4e/cuTp06hUmTJnn7kIIOP8yPGloEgjpvBw37ZMdP2rhxIzZu3Ahg4MaT/rZNSEgYcJ+OqqqqXCxhaKLLZ6GhhhZCPPI4KdKIFmGgpEiIR/h7ijT2WSAslBQJ8QTfT9FCSVEY6J4iIR6hSWaFhlqfCfEINbQIDd1TJMQj1E9RaCgpEuIR+5mi2cpgtnq3BZqSoj9QQwshHglzGMPu7bNFSor+QA0thHhELhXDPkzc2/cVKSn6AzW0EOIRkUj0eKJZE10+Bz+6pxg0HBeu4jgOSUlJ+Pjjj2GxWFBVVQWRSMRPD+Zo/Pjx2L1794CPied81VeRxj77Az9LDl0+B4P58+fj4MGDMBqNOHXqFHQ6HWQyGbRarb+LFtIUPpoUgs4U/YGWIwgq9oWrEhISkJeXh4yMDPz444/+LlbIs3fg9vYyp3Sm6A/U0AIwBpi7hr+dzda7nUkCuDHzNmRKwMOJZhUKRZ91msnI4/sq0uWzANjXaAnlM0VzF/BJ3LA3EwMY7cn7bmoCuKfc2pQxhsrKSpSXl2P16tX8888880yf2K4uNxI+GZbHDS2UFIMfLUcQVOyr+ZnNZthsNixduhTbt2/H77//DgA4f/58n4Wk0tLS/FDS0BLmo6F+lBT9gW99DuHLZ5my96xtmGw2G9o7OhCucm/hKsiUw97Evpofx3GIi4uDVOr8tUlMTMTo0c7nr0/GEO/z1fhn+sv5g4X6KUIkcu8y1mYDZNbebd1Jim5wXM2PBA5frf1Mrc/+QP0UCfFYmLQ3KRot1Hk7+FHrMyEe89WZIl0++wMN8wsagy0ylZaWNuCiVTdv3hz0MfGcrxpa6EzRH+jymRCP+aqhhZLiSLPZAJul9/dQbn0mxEMKrjd9ebufIiXFkWY/SwTo8pkQD9CZolA4JkVqaCHEbXRPUSjso1kAQExnioS4i/opCoV93LNYOmKdjwOFzebd/mSkr1CqY3s/RW8vR0BdckZaCA7x4zgOYrEYTU1NGDduHDiOg8jNmWpsNhtMJhN6enrcG+YnUIwxmEwm3LlzB2KxGBwn/J4Nj1f08+4/AkqKIy0E51IUi8VITExEc3MzmpqGP97ZEWMM3d3dUCgUbidWIVMqldBoNCHxD4MmhBAKftZt4f8nd8RxHDQaDSwWC6xW9z/EZrMZNTU1mDNnDmSy0PnH4gqJRAKpVBoy/yxoQgihCOEhfiKRCDKZzKNkJpFIYLFYEBYWRkkxxPGXz0JqaCkuLsbMmTOhUqkQHR2NnJwcNDQ0OMX09PRAp9NhzJgxGDVqFF5//XW0tLQ4xdy+fRvZ2dlQKpWIjo7GunXrYLFYnGKqqqowffp0yOVyJCUl9Tt8a+/evRg/fjzCwsKQmpqK3377zevHHIqXz4T4giD7KVZXV0On0+HixYuoqKiA2WxGZmYmOjs7+Zi1a9fip59+wtGjR1FdXY2mpiYsWrSIf91qtSI7OxsmkwkXLlzA4cOHcejQIWzdupWPuXHjBrKzs5Geno66ujoUFBRg5cqVKC8v52PKyspQWFiIbdu24fLly0hOTkZWVhZaW1u9e9DW0Lx8JsTb7EnRYmMwW73Y2MICSGtrKwPAqqurGWOMPXz4kMlkMnb06FE+pr6+ngFger2eMcbYqVOnmFgsZgaDgY/Zv38/Cw8PZ0ajkTHG2Pr169mLL77o9F6LFy9mWVlZ/ONZs2YxnU7HP7ZarSwuLo4VFxe7VPa2tjYGgB05coSZTKaBA69WMLYtnLH9/+PSfokzk8nEjh8/PngdE48ESx33mC0sYcNJlrDhJGvrHrqs9u9oW1vboHEBdU+xra0NABAVFQUAuHTpEsxmMzIyMviY559/HhqNBnq9HrNnz4Zer8dLL72EmJgYPiYrKwt5eXn466+/8PLLL0Ov1zvtwx5TUFAAADCZTLh06RKKior418ViMTIyMqDX6/stq9FohNFo7FP2rq4u3Lt3b8D7XaJ79yA1MtiMIlhp8aNhM5vNQ9Yx8Uyw1DFjDDB1wcaAJsMdmFWD36fv6Oh4vN0gAiYp2mw2FBQU4NVXX8WUKVMAAAaDARzH9ZnqPSYmBgaDgY9xTIj21+2vDRbT3t6O7u5uPHjwAFartd+Yf/75p9/yFhcXY8eOHX2eX7lypYtHfB4oGOtiLCFkMJN3ux7b0dGBiIiIAV8PmKSo0+lw5coV/PLLL/4uikuKiopQWFjIP7bZbLh16xamTZuGxsZGhIeH+7F0wtXe3o74+HiqYx8Sah0zxtDR0YG4uMFXkQyIpJifn4+TJ0+ipqbGablItVoNk8mEhw8fOp0ttrS0QK1W8zFPthLbW6cdY55ssW5paUF4eDgUCgUkEgkkEkm/MfZ9PEkul0Mudz5dt3eYDQ8PF9SHKRBRHfueEOt4sDNEO7+2PjPGkJ+fj2PHjuHs2bNITEx0ej0lJQUymQyVlZX8cw0NDbh9+za0Wi0AQKvV4s8//3RqJa6oqEB4eDheeOEFPsZxH/YY+z44jkNKSopTjM1mQ2VlJR9DCAkR3mgFcldeXh6LiIhgVVVVrLm5mf/p6uriY95//32m0WjY2bNnWW1tLdNqtUyr1fKvWywWNmXKFJaZmcnq6urY6dOn2bhx41hRUREfc/36daZUKtm6detYfX0927t3L5NIJOz06dN8TGlpKZPL5ezQoUPs77//ZqtWrWKjR492atUeiqutW8R9VMe+F+p17NekCKDfn4MHD/Ix3d3d7IMPPmCRkZFMqVSyhQsXsubmZqf93Lx5ky1YsIApFAo2duxY9tFHHzGz2ewUc+7cOTZt2jTGcRybMGGC03vYffHFF0yj0TCO49isWbPYxYsXh3U8PT09bNu2baynp2dY2xHXUR37XqjXsYixIdqnCSEkhAh/Kg1CCBkGSoqEEOKAkiIhhDigpEgIIQ4oKXrJiEw7FsJqamrw2muvIS4uDiKRCMePH/d3kQTHlan8QgElRS8YsWnHQlhnZyeSk5Oxd+9efxdFsFyZyi8UUJccL0hNTcXMmTOxZ88eAL2jYeLj47F69Wps3LjRz6UTHpFIhGPHjiEnJ8ffRRG0O3fuIDo6GtXV1ZgzZ46/izNi6EzRQ/ZpxxynJhtq2jFCgsGTU/mFCkqKHrp79+6A047Zpy4jJNj0N5VfqAiIWXIIIYEl2Kby8yZKih4aO3bssKcdIySQDTSVX6igy2cP0bRjRCjYEFP5hQo6U/SCwsJCrFixAjNmzMCsWbOwe/dudHZ24q233vJ30QTj0aNHuHbtGv/4xo0bqKurQ1RUFDQajR9LJhw6nQ5HjhzBiRMnoFKp+HviERERUCgUfi7dCPLrHD0C4um0Y2Rw586d63eauRUrVvi7aILRX/3iian8QgH1UySEEAd0T5EQQhxQUiSEEAeUFAkhxAElRUIIcUBJkRBCHFBSJIQQB5QUCSHEASVFQghxQEmREEIcUFIkhBAHlBQJIcQBJUVCCHHw/wmR/S1k2M2WAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 300x600 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "jkivwCHtA1lt",
        "colab_type": "text"
      },
      "source": [
        "#Calculating Pore Pressure Using d-exponent"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "WyIgN-NI8Ec6",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 35
        },
        "outputId": "eacb45f3-f00e-481f-c0a7-195eb2a61422"
      },
      "source": [
        "d1 = np.array([8100,9000,9600,10100,10400,10700,10900,11100,11300,11500,11600,11800,12100,12200,12300,\n",
        "               12700,12900,13000,13200,13400,13500,13600,13700,13800,13900,14000,\n",
        "               14200,14400,14600,14800,14900,15000,15200,15300,15400,15500,15700,16200,16800])\n",
        "\n",
        "dc = np.array([1.52,1.55,1.57,1.49,1.58,1.60,1.61,1.57,1.64,1.48,1.61,1.54,1.58,1.67,1.41,1.27,1.18,1.13,1.22,\n",
        "               1.12,1.12,1.07,1.00,0.98,1.00,0.91,0.93,0.86,0.80,0.86,0.80,0.90,0.82,0.87,0.92,\n",
        "               0.87,0.80,0.80,0.65])\n",
        "\n"
      ],
      "execution_count": 17,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "39"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 17
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "YP5V90zrExTt",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "y = d1[d1<=12500]\n",
        "x = dc[0:15]\n",
        "\n",
        "# z = np.polyfit(x, y, 1)\n",
        "# p = np.poly1d(z)\n",
        "d11 = np.array([6000,8100,9000,9600,10100,10400,10700,10900,11100,11300,11500,11600,11800,12100,12200,12300,\n",
        "               12700,12900,13000,13200,13400,13500,13600,13700,13800,13900,14000,\n",
        "               14200,14400,14600,14800,14900,15000,15200,15300,15400,15500,15700,16200,16800])\n",
        "dmod = 1.15 + 0.000038*d11"
      ],
      "execution_count": 58,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "EXv6k2-yFJmT",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 35
        },
        "outputId": "67248fd5-cc03-416a-b5c2-d5ab34c6b187"
      },
      "source": [
        "#dmod = 1.15 + 0.000038D"
      ],
      "execution_count": 39,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "12"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 39
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "4hN3QXp7CpfQ",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 545
        },
        "outputId": "167879fe-8ae3-4ddb-b8d9-e60897dd3482"
      },
      "source": [
        "plt.style.use('default')\n",
        "plt.figure(figsize=(3,6))\n",
        "\n",
        "plt.scatter(dc,d1,label='d-exp(normal)')\n",
        "plt.scatter(x,y,label='d(normal)-linear trend')\n",
        "plt.plot(dmod,d11,'red',label='d-exp(corrected)')\n",
        "plt.ylim(16800,5000)\n",
        "\n",
        "plt.legend(loc='upper left')\n",
        "plt.title('Pore Pr. Calc. with D-Exponent')\n",
        "plt.ylabel('Depth(ft)')\n",
        "plt.grid()"
      ],
      "execution_count": 67,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAU8AAAIQCAYAAAAB/exNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeVxU1fsH8M+wDPui7CgCLoAporkQKqKJQpql5U6KhmnlkrmlmQtpbmhf03LLX2rlmmZZbuCWpESGC4sCirizqOzIPs/vj2luXEEZcGAGeN6v17x0zj333HPvzDzc9TwSIiIwxhirFi11d4AxxuojDp6MMVYDHDwZY6wGOHgyxlgNcPBkjLEa4ODJGGM1wMGTMcZqgIMnY4zVAAdPxhirAQ6eTKNs374dEokEt27dUndXRBYvXgyJRFKtuo8eParlXjF14uBZDYoftuKlr68PFxcXTJkyBWlpaerunqB3796ifjZt2hRdu3bFd999B5lMptJlFRYW4n//+x88PT1hZmYm2iaJiYkqXZamWbZsGX755ReVtztu3DjR52dsbIyWLVti6NChOHDgQLU+QycnJ1Fb5V/+/v4q73t9s2HDBmzfvr1G8+qotiuNw+effw5nZ2cUFhbizz//xMaNG3HkyBHExsbC0NBQ3d0DADRv3hzLly8HADx8+BDff/89goKCkJiYiBUrVqhkGY8ePYK/vz+ioqLw+uuvY/To0TA2NkZCQgL27NmDLVu2oLi4WCXLUrfPPvsMc+fOFZUtW7YMQ4cOxeDBg1W+PD09PWzduhUAUFBQgNu3b+O3337D0KFD0bt3b/z6668wNTVVqq2OHTti5syZFcrt7e1V2uf6aMOGDbC0tMS4ceOqPzMxpW3bto0A0IULF0TlM2bMIAC0a9euF2pfJpPRkydPXqgNIiIfHx9q166dqCw/P5+aN29ORkZGVFxcXOl8ZWVlVFBQoPRyBg4cSFpaWrR///4K0woLC2nmzJnV6zj9t42Tk5OrPW9dMzIyosDAwArlixYtIgD08OHDGrUbGBhIRkZGlU5bvnw5AaDhw4cr1ZajoyMNHDiwRv1oDNq1a0c+Pj41mpcP21Xg1VdfBQAkJycDAEpLS7FkyRK0atUKenp6cHJywqeffoqioiLRfE5OTnj99ddx/PhxdOnSBQYGBti8eTMAICsrC9OnT4eDgwP09PTQunVrrFy5ssaH3YaGhnjllVeQn5+Phw8fAgAkEgmmTJmCnTt3ol27dtDT08OxY8eUai8yMhKHDx9GUFAQ3n777QrT9fT0sHr1auF9dHQ0xo0bh5YtW0JfXx+2trZ499138fjxY6WWd/ToUfj4+MDExASmpqbo2rUrdu3apdS8CkQES0tLzJgxQyiTyWQwNzeHtrY2srKyhPKVK1dCR0cHeXl5ACqe85RIJMjPz8eOHTuEw+Cn916ysrIwbtw4mJubw8zMDOPHj8eTJ0+q1eenzZ07F/3798dPP/2kstMi6enpsLKyQu/evUHlBlm7ceMGjIyMMGLECKGsd+/eaN++PaKiotC9e3cYGBjA2dkZmzZtqrTdoKAg2NjYQF9fHx4eHtixY4eozq1btyCRSLB69Wps2bJF+M107doVFy5cqNBmfHw8hg4diqZNm0JfXx9dunTBoUOHRHUUp9fOnTuHGTNmwMrKCkZGRhgyZIjw3Qfkv7+4uDj88ccfwmfYu3dvpbcbH7arQFJSEgDAwsICADBhwgTs2LEDQ4cOxcyZMxEZGYnly5fj2rVrOHjwoGjehIQEjBo1CpMmTcJ7770HV1dXPHnyBD4+Prh//z4mTZqEFi1a4Pz585g3bx5SUlKwdu3aGvXz5s2b0NbWhrm5uVB26tQp7Nu3D1OmTIGlpSWcnJyUakvxhR0zZoxS9cPCwnDz5k2MHz8etra2iIuLw5YtWxAXF4e//vrruRdjtm/fjnfffRft2rXDvHnzYG5ujkuXLuHYsWMYPXq0UssH5AGvR48eOHv2rFAWHR2N7OxsaGlp4dy5cxg4cCAAIDw8HJ06dYKxsXGlbf3www+YMGECunXrhokTJwIAWrVqJaozfPhwODs7Y/ny5bh48SK2bt0Ka2trrFy5Uuk+V2bMmDEIDQ1FWFgYXFxcqqxfUlJS6cUrIyMjGBgYwNraGhs3bsSwYcOwfv16TJs2DTKZDOPGjYOJiQk2bNggmi8zMxMDBgzA8OHDMWrUKOzbtw8ffPABpFIp3n33XQDyUw29e/fGjRs3MGXKFDg7O+Onn37CuHHjkJWVhY8++kjU5q5du5Cbm4tJkyZBIpFg1apVeOutt3Dz5k3o6uoCAOLi4tCjRw80a9YMc+fOhZGREfbt24fBgwfjwIEDGDJkiKjNqVOnokmTJli0aBFu3bqFtWvXYsqUKdi7dy8AYO3atZg6dSqMjY0xf/58AICNjY2SnwL4sL06FIeUJ06coIcPH9Ldu3dpz549ZGFhQQYGBnTv3j26fPkyAaAJEyaI5p01axYBoFOnTglljo6OBICOHTsmqrtkyRIyMjKixMREUfncuXNJW1ub7ty589x++vj4kJubGz18+JAePnxI165do2nTphEAGjRokFAPAGlpaVFcXFy1t8WQIUMIAGVmZipVv7LTEbt37yYAdPbsWaHs6cP2rKwsMjExIU9PzwqnFGQyWbX7HRISQtra2pSTk0NEROvWrSNHR0fq1q0bffLJJ0QkP31hbm5OH3/8sTCf4lC8vKoO2999911R+ZAhQ8jCwqLKPj7vsJ2I6NKlSwRA1L9nUXzHKnstX75cVHfUqFFkaGhIiYmJFBISQgDol19+EdXx8fEhALRmzRqhrKioiDp27EjW1tbCKaG1a9cSAPrxxx+FesXFxeTl5UXGxsbC9k9OTiYAZGFhQRkZGULdX3/9lQDQb7/9JpT17duX3N3dqbCwUCiTyWTUvXt3atOmjVCm+A75+vqKviMff/wxaWtrU1ZWllDGh+11zNfXF1ZWVnBwcMDIkSNhbGyMgwcPolmzZjhy5AgAiA4NAQgn7A8fPiwqd3Z2hp+fn6jsp59+gre3N5o0aYJHjx4JL19fX5SVlYn2nJ4lPj4eVlZWsLKyQtu2bbF+/XoMHDgQ3333naiej48PXnrppWpvg5ycHACAiYmJUvUNDAyE/xcWFuLRo0d45ZVXAAAXL1585nxhYWHIzc3F3Llzoa+vL5qm7K1D5Xl7e6OsrAznz58HIN/D9Pb2hre3N8LDwwEAsbGxyMrKgre3d7XbL+/999+vsOzHjx8L266mFHvDubm5StX39PREWFhYhdeoUaNE9b7++muYmZlh6NChWLBgAcaMGYM333yzQns6OjqYNGmS8F4qlWLSpElIT09HVFQUAODIkSOwtbUVLUNXVxfTpk1DXl4e/vjjD1GbI0aMQJMmTYT3im1/8+ZNAEBGRgZOnTqF4cOHIzc3V/hNPH78GH5+frh+/Tru378vanPixImi74jis799+7ZS260qfNheA9988w1cXFygo6MDGxsbuLq6QktL/nfo9u3b0NLSQuvWrUXz2NrawtzcvMIH5+zsXKH969evIzo6GlZWVpUuPz09vco+Ojk54dtvvxVuqWrTpg2sra0r1Kts+cpQXOnNzc0VnQZ4loyMDAQHB2PPnj0V+p+dnf3M+RSnRNq3b1+jfj7t5ZdfhqGhIcLDw+Hn54fw8HAEBwfD1tYW69evR2FhoRBEe/bs+ULLatGihei9IjhkZmYqfaW8MorzsIo/XNnZ2SgoKBCmS6VSNG3aVHhvaWkJX1/fKttt2rQp1q1bh2HDhsHGxgbr1q2rtJ69vT2MjIxEZYrTB7du3cIrr7yC27dvo02bNsLvQqFt27YAUOF38LxtBcjPvxIRFixYgAULFlTar/T0dDRr1kzpNl8UB88a6NatG7p06fLcOsruFZXfI1OQyWTo168f5syZU+k8ypznMjIyUuoHU9nyleHm5gYAiImJUWoPbfjw4Th//jxmz56Njh07wtjYGDKZDP7+/iq/9/R5dHV14enpibNnz+LGjRtITU2Ft7c3bGxsUFJSgsjISISHh8PNze2Zf7yUpa2tXWk5vWDmm9jYWAAQ/kB/9NFHogsxPj4+OHPmTI3aPn78OAB5gLl3755SfxhVoaptpfiOzJo1q8KRmsLTOyy1tf0VOHiqmKOjI2QyGa5fvy78lQWAtLQ0ZGVlwdHRsco2WrVqhby8PKWCn7oMGjQIy5cvx48//lhl8MzMzMTJkycRHByMhQsXCuXXr1+vcjmKizCxsbEVfhw15e3tjZUrV+LEiROwtLSEm5sbJBIJ2rVrh/DwcISHh+P111+vsp2anDZQhR9++AESiQT9+vUDAMyZMwfvvPOOML384W91HDt2DFu3bsWcOXOwc+dOBAYGIjIyEjo64jDx4MED5Ofni/Y+FVf+FRccHR0dER0dDZlMJtr7jI+PF6ZXR8uWLQHI//ip8nfxIp8hn/NUsQEDBgBAhSviX375JQAIV3OfZ/jw4YiIiBD2AsrLyspCaWmpCnpauezsbMTHxz/3UBoAvLy84O/vj61bt1b6lE1xcTFmzZoF4L89gKf/4itz10D//v1hYmKC5cuXo7CwUDStfHuPHj1CfHy8UrcCeXt7o6ioCGvXrkXPnj2FH5C3tzd++OEHPHjwQKm9aSMjI9HtTXVhxYoVCA0NxYgRI9CmTRsAwEsvvQRfX1/h1blz52q3m5WVJdw9sGzZMmzduhUXL17EsmXLKtQtLS0VbqkD5J/15s2bYWVlJSx7wIABSE1NFa5sK+Zbv349jI2N4ePjU63+WVtbo3fv3ti8eTNSUlIqTC9/C1J1vMhnyHueKubh4YHAwEBs2bIFWVlZ8PHxwd9//40dO3Zg8ODB6NOnT5VtzJ49G4cOHcLrr7+OcePGoXPnzsjPz0dMTAz279+PW7duwdLSslb6f/DgQYwfPx7btm2r8qmL77//Hv3798dbb72FQYMGoW/fvjAyMsL169exZ88epKSkYPXq1TA1NUWvXr2watUqlJSUoFmzZggNDRXui30eU1NT/O9//8OECRPQtWtXjB49Gk2aNMGVK1fw5MkT4XD166+/RnBwME6fPl3lvXpeXl7Q0dFBQkKCcJsRAPTq1QsbN24EAKWCZ+fOnXHixAl8+eWXsLe3h7OzMzw9PaucTxmlpaX48ccfAcgvsN2+fRuHDh1CdHQ0+vTpgy1btijd1v3794W2yjM2Nhaejvroo4/w+PFjnDhxAtra2vD398eECROwdOlSvPnmm/Dw8BDms7e3x8qVK3Hr1i24uLhg7969uHz5MrZs2SLcVjRx4kRs3rwZ48aNQ1RUFJycnLB//36cO3cOa9euVfpCY3nffPMNevbsCXd3d7z33nto2bIl0tLSEBERgXv37uHKlSvVbrNz587YuHEjli5ditatW8Pa2lq4b7tKNbpG30g96wmjp5WUlFBwcDA5OzuTrq4uOTg40Lx580S3WBA9/+mP3NxcmjdvHrVu3ZqkUilZWlpS9+7dafXq1c98QkihsieMKgOAJk+eXOk6btu2rcr5ieS3IK1evZq6du1KxsbGJJVKqU2bNjR16lS6ceOGUO/evXs0ZMgQMjc3JzMzMxo2bBg9ePCAANCiRYsqLP/pJ4wOHTpE3bt3JwMDAzI1NaVu3brR7t27hemK24NOnz6tVL+7du1KACgyMlLURwDk4OBQoX5ltyrFx8dTr169yMDAgAAIty096wkjZZ+eCgwMFN1SZGhoSE5OTvT222/T/v37qaysTKl1JHr+rUqOjo5E9N9tQeVvPyIiysnJIUdHR/Lw8BC+c4rv1j///ENeXl6kr69Pjo6O9PXXX1dYdlpaGo0fP54sLS1JKpWSu7t7he+V4lalkJCQCvM//d0gIkpKSqKxY8eSra0t6erqUrNmzej1118XPeX2rN/p6dOnK3xHUlNTaeDAgWRiYkIAqnXbkuTfTjLGWJV69+6NR48eCRetGjM+58kYYzXAwZMxxmqAgydjjNUAn/NkjLEa4D1PxhirAQ6ejDFWA3yTvIrIZDI8ePAAJiYmantsjzH2bESE3Nxc2NvbVxiwpCY4eKrIgwcP4ODgoO5uMMaqcPfuXTRv3vyF2+HgqSKKx83u3r37QsONNRYlJSUIDQ1F//79hUf6WO1q7Ns8JycHDg4ONXo0tDIcPFVEcahuamrKwVMJJSUlMDQ0hKmpaaP8IasDb3M5VZ1W4wtGjDFWAxw8GWOsBjh4MsZYDfA5zzpWVlaGkpISdXdD7UpKSqCjo4PCwkKUlZWpuzuNQn3Z5lKpVCW3EtU2Dp51hIiQmppa5yOPayoigq2tLe7evcv3xdaR+rLNtbS04OzsDKlUqu6uPBcHzzqiCJzW1tYwNDTU6C9vXZDJZMjLy4OxsXG92MtoCOrDNlc8bJKSkoIWLVpo9O+Eg2cdKCsrEwKnhYWFurujEWQyGYqLi6Gvr6+xP+SGpr5scysrKzx48AClpaUafUuV5m7BBkRxjtPQ0FDNPWFM8ykO1zX5vCzAwbNOafIhCGOaor78Tjh4MsZYDXDwZNXWu3dvTJ8+Xd3dQHFxMVq3bo3z58+ruytq5+TkhLVr1wKQbxcnJyf8888/au5Vw8bBk9VbmzZtgrOzM7p3767urtRI+YCnSlKpFLNmzcInn3yi8rbZfzh4snqJiPD1118jKCioVpdT2QMNxcXFtbpMVQgICMCff/6JuLg4dXelweLgyZ4rPz8fY8eOhbGxMezs7LBmzZoq58nKysKECRNgZWUFU1NTvPrqq7hy5QoA4OHDh7C1tcXy5cuF+ufPn4dUKsXJkycBAIsXL0bHjh2xefNmODg4wNDQEMOHD0d2drYwT1RUFJKSkjBw4EDRsu/du4dRo0ahadOmMDIyQpcuXRAZGSlM37hxI1q1agWpVApXV1f88MMPovklEgk2btyIN954A0ZGRvjiiy+E/mzduhXOzs7Q19evcj0VfvvtN3Tt2hX6+vqwtLTEkCFDAMhPfdy+fRsff/wxJBKJ6CLJn3/+CW9vbxgYGMDBwQHTpk1Dfn6+MD09PR2DBg2CgYEBnJ2dsXPnzgqfQZMmTdCjRw/s2bOnys+L1QwHT3UhAvLz1fOqRs6/2bNn448//sCvv/6K0NBQnDlzBhcvXnzuPMOGDUN6ejqOHj2KqKgovPzyy+jbty8yMjJgZWWF7777DsHBwbh06RJyc3MxZswYTJkyBX379hXauHHjBvbt24fffvsNx44dw6VLl/Dhhx8K08PDw+Hi4iIamzEvLw8+Pj64f/8+Dh06hCtXrmDOnDmQyWQAgIMHD+Kjjz7CzJkzERsbi0mTJmH8+PE4ffq0qP+LFy/GkCFDEBMTg3fffVfoz4EDB/Dzzz/j8uXLVa4nABw+fBhDhgzBgAEDcOnSJZw8eRLdunUDAPz8889o3rw5Pv/8c6SkpCAlJQUAkJSUBH9/f7z99tuIjo7G3r178eeff2LKlClC/8aNG4e7d+/i9OnT2L9/PzZs2ID09PQKn0O3bt0QHh5exSfMaoyYSmRnZxMAys7OrjCtoKCArl69SgUFBf8V5uURycNY3b/y8pRap9zcXJJKpbRv3z6h7PHjx2RgYEAfffRRpfOEh4eTqakpFRYWispbtWpFmzdvFt5/8MEH1Lp1axo1ahS5u7uL6i9atIi0tbXp3r17QtnRo0dJS0uLUlJSiIjoo48+oldffVW0jM2bN5OJiQk9fvy40r51796d3nvvPVHZsGHDaMCAAcJ7ADR9+nRRnUWLFpGuri6lp6dXaz29vLwoICCg0r4QETk6OtL//vc/UVlQUBBNnDhRVBYeHk5aWlpUUFBACQkJBID+/vtvYfq1a9cIQIW2vvrqK3JychLel5WVUWZmJpWVlT2zT5qg0t+LCjzvN1oTvOfJnikpKQnFxcXw9PQUypo2bQpXV1cAwLJly2BsbCy87ty5gytXriAvLw8WFhaiacnJyUhKShLaCQkJQWlpKfbv34+dO3dCT09PtOwWLVqgWbNmwnsvLy/IZDIkJCQAAAoKCoTDZ4XLly+jU6dOaNq0aaXrc+3aNfTo0UNU1qNHD1y7dk1U1qVLlwrzOjo6wsrKSnivzHpevnxZtDetjCtXrmD79u2iNv38/CCTyZCcnIxr165BR0cHnTt3FuZxc3ODubl5hbYMDAzw5MmTai2fKY8fz1QXQ0MgL099y1aB999/H8OHDxfe29vbIy8vD3Z2djhz5kyF+uV/4ElJSUhNTYVMJsOtW7fg7u5erWVbWloiJiZGVGZgYFC9FXgGIyOjKsuUWc+a9CcvLw+TJk3CtGnTKkxr0aIFEhMTlW5LcZqE1Q4OnuoikQCV/Eg1SatWraCrq4vIyEi0aNECAJCZmYnExET4+PigadOmFfbyXn75ZaSmpkJHRwdOTk6VtltcXIyxY8diyJAhaN++PSZMmICYmBhYW1sLde7cuYMHDx7A3t4eAPDXX39BS0tL2Ovt1KkTNm7cCCISLrZ06NABW7duRUZGRqV7n23btsW5c+cQGBgolJ07dw4vvfRStbeNMuvZoUMHnDx5EuPHj690ulQqrfAI4ssvv4yrV6+idevWlc7j5uaG0tJSREVFoWvXrgCAhISESkfrio2NRadOnaqxVqxaVHLwz6p/zrOeeP/998nR0ZFOnjxJMTEx9MYbb5CxsfEzz3nKZDLq2bMneXh40PHjxyk5OZnOnTtHn376KV24cIGIiGbNmkVOTk50+/ZtKikpoZ49e9LAgQOFNhYtWkRGRkbk6+tLly9fprNnz5KLiwuNHDlSqPPo0SPS1dWlmJgYoayoqIhcXFzI29ub/vzzT0pKSqL9+/fT+fPniYjo4MGDpKurSxs2bKDExERas2YNaWtr0+nTp4U2ANDBgwdF67Ro0SLy8PCo9nqePn2atLS0aOHChXT16lWKjo6mFStWCG3069eP3njjDbp37x49fPiQiIiuXLlCBgYGNHnyZLp06RIlJibSL7/8QpMnTxbm8/f3p06dOtFff/1F//zzD/Xs2ZMMDAwqnPN0dHSk77//XnjP5zxVe86Tg6eKNNTgmZubS++88w4ZGhqSjY0NrVq1inx8fJ4ZPImIcnJyaOrUqWRvb0+6urrk4OBAAQEBdOfOHTp9+jTp6OjQH3/8IfyQk5OTydTUlDZs2EBE/wWrDRs2kL29Penr69PQoUMpIyNDtJzhw4fT3LlzRWW3bt2it99+m0xNTcnQ0JC6dOlCkZGRwvQNGzZQy5YtSVdXl1xcXETBhUj54FnVeiocOHCAOnbsSFKplCwtLemtt94SpkVERFCHDh1IT0+Pyu/H/P3339SvXz8yNjYmIyMj6tChA33xxRfC9JSUFBo4cCDp6elRixYt6Pvvv69w8en8+fNkbm5OT548Eco4eKo2eEqIqnHfCnumnJwcmJmZITs7u0L2zMLCQiQnJ4vuEWzsZDIZcnJyYGpqWmF4tMWLF+OXX34Rbgl6lujoaPTr1w9JSUkwNjauze7WOyNGjICHhwc+/fRToex521yT1Nbv5Xm/0ZrQ3C3IWBU6dOiAlStXIjk5Wd1d0SjFxcVwd3fHxx9/rO6uNGh8wYjVa+PGjVN3FzSOVCrFZ599pu5uNHi858k0zuLFi6s8ZGdM3Th4MsZYDXDwZIyxGuDgyRhjNcDBkzHGaoCDJ2OM1QAHT8YYqwEOnqzaOAGc6pw5cwYSiUQY2OPYsWPo2LGjMIAz01waHzzv37+Pd955BxYWFjAwMIC7u7soKyARYeHChbCzs4OBgQF8fX1x/fp1URsZGRkICAiAqakpzM3NERQUhLynhoOLjo6Gt7c39PX14eDggFWrVtXJ+rGaq+8J4Crj7+8PXV3dSlNrMM2i0cEzMzMTPXr0gK6uLo4ePYqrV69izZo1aNKkiVBn1apVWLduHTZt2oTIyEgYGRnBz88PhYWFQp2AgADExcUhLCwMv//+O86ePYuJEycK03NyctC/f384OjoiKioKISEhWLx4MbZs2VKn66uMMhkhIukxfr18HxFJj1Ema5xDE1AdJYAD6j7h27hx47Bu3bo6XSarAZUML1JLPvnkE+rZs+czp8tkMrK1taWQkBChLCsri/T09Gj37t1ERHT16lUCIAwTRiRP6SCRSOj+/ftEJB9pp0mTJlRUVCRatqurq9J9rYtRlY7GPKBXlp0gx09+F16vLDtBR2MevFC7z5OXl0djxowhIyMjsrW1pdWrV1c5qlJmZiYFBQWRpaUlmZiYUJ8+fejy5ctERJSenk42Nja0dOlSYYSfc+fOka6uLp04cYKI/hvFaNOmTdS8eXMyMDCgYcOGUVZWlrCMCxcukJaWFuXk5AhlycnJBIAOHDhAvXv3JgMDA+rQoYMwJJ3C/v376aWXXiKpVEqOjo60evVq0XRHR0f6/PPPacyYMWRiYkKBgYG0bds2MjMzo99++41cXFzIwMCA3n77bcrPz6ft27eTo6MjmZub09SpU6m0tFRo6/vvv6fOnTuTsbEx2djY0KhRoygtLU2Yfvr0aQJAmZmZQtnt27cJAN24caM6H1WVeFSlRpSG49ChQ+jSpQuGDRsGa2trdOrUCd9++60wPTk5GampqfD19RXKzMzM4OnpiYiICABAREQEzM3NRakVfH19oaWlJWRVjIiIQK9evSCVSoU6fn5+SEhIQGZmZqV9KyoqQk5OjugFyFPVVvYiIshkshq/jkQ/wAc/XkRKdqGoH6nZhfjgx4s4Ev3ghdp/1mvWrFn4448/cPDgQRw7dgynT5/GxYsXn7s+Q4cORVpaGg4fPowLFy6gU6dO6Nu3Lx49egQLCwts3bpVSACXk5ODMWPGYPLkyejTpw9kMhmISEgA9+uvv+LIkSO4dOkSPvjgA2EZZ8+ehYuLC4yMjETLBoD58+djxowZuHjxItq0aYNRo0ahuLgYMpkMFy5cwPDhwzFixAhcuXIFCxcuxIIFC/Ddd9+J2li9ejU6dOiAqKgozJ8/HzKZDE+ePMFXX32FXbt24ciRIzhz5gwGDx6Mw4cP4/fff8eOHTuwefNm7Nu3T2irqKhIWNeff/4ZtzeFKW0AACAASURBVG7dQmBgYIU+y2QyyEpKQHfvwqGkBDY2Nvjjjz9U+lnSvwOoveh3sS5eRPTM39KLvFRJowcGuXnzJjZu3IgZM2bg008/xYULFzBt2jRIpVIEBgYiNTUVAGBjYyOaz8bGRpiWmpoqGqEcAHR0dNC0aVNRHWdn5wptKKaVP02gsHz5cgQHB1coDw0NheFTaS50dHRga2uLvLy8Gh0ClskIwb/FobIDdAIgARD8Wxw8mxtAW0tSSa2aycvLw3fffYfNmzcLo5avX78e7dq1Q3FxsfAHo7yIiAj8/fffuH79upCXaMGCBTh48CB+/PFHjBs3Dj179sTYsWMxceJEdOzYEfr6+pg7d67QXlFREQoLC7F+/XphJPnly5djxIgRWLRoEWxsbHD9+nVYW1uL+qA4j/3hhx/C29sbADBr1ix4eXnh8uXLcHFxwapVq+Dj4yOkuXjrrbdw+fJlhISE4K233gIgD2Te3t6YMGGC0HZhYSFKSkqwatUq4bsyaNAg7Nu3DwkJCTA2Nkbz5s3h7e2N0NBQvPbaawCAoUOHCm1YWlriiy++wKuvvooHDx7A2NhYyDGUm5sLLQBmaWmQALCxtsb169cr3cYvKjc3V+VtqlJxcTEKCgpw9uxZlJaWqqxdVedz0ujgKZPJ0KVLFyxbtgyAPPVCbGwsNm3aJEqloA7z5s3DjBkzhPc5OTlwcHBA//79Kx3P8+7duzA2Nq7R+IR/3XyMtNxnB10CkJZbjISMUrzS0qLa7T9LcnIyiouL0bt3b2GdTE1N4erqCqlUim+++UaUfz02NhZJSUnIz89Hq1atRG0VFBTgwYMHQjtr166Fu7s7fv31V1y4cEGUa0dPTw8tWrSAm5ubUNa3b1/IZDI8ePAAbdq0QVlZGYyMjETbWjGmZ7du3YTyNm3aAJD/cExNTZGUlIQ33nhDNF+fPn2wadMmGBkZQVtbG1paWnjllVdEdfT19WFoaAgPDw+hzMHBAU5OTkKAB+R5nLKysoR5o6KiEBwcjOjoaGRmZgp7mllZWbC3txf+0JqYmMDU3BzQ0wOKimBsaIiysjKVjDupQETIzc2FiYmJKE+8piksLISBgQF69eql8vE8VUmjg6ednV2F/DJt27bFgQMHAAC2trYAgLS0NNjZ2Ql10tLS0LFjR6HO0zmtS0tLkZGRIcxva2uLtLQ0UR3Fe0Wdp+np6VXI+AgAurq60NXVFZWVlZVBIpFAS0urRoPQPsxTbm/1YV6xSge5VbRVWb8lEgk++OADjBgxQihr3rw58vPzn5sYTdHOzZs3hQRwd+7cEQUlxQ+7/DKf7ouVlRViY2MrraOnpyf8X1tbu8J0xWfxvPU0NjauUEdXV1epMiKClpYW8vPz8dprr8HPzw87d+6ElZUV7ty5Az8/P5SWloqWJ/zfwAAoKkJGRgasra1V+nkqAvfT669ptLS0IJFIKv0tvQhVtgVo+NX2Hj16CKlmFRITE+Ho6AgAcHZ2hq2tLU6ePClMz8nJQWRkJLy8vADIU9ZmZWUhKipKqHPq1CnIZDIhpa6XlxfOnj0rOicSFhYGV1fXSg/Z65q1iXJ/fZWtp6zyCeAUFAngAHka4tatWwsvHR0dUWK08tNat24NS0tLAOIEcJ9//jkmTJhQ4Q+cIgGcQmUJ4OLj44XzeMpSJIEr79y5c3BxcREFWlWIj4/H48ePsWLFCnh7e8PNza3CelZgaIjCoiIk3brFyds0nEYHz48//hh//fUXli1bhhs3bmDXrl3YsmULJk+eDED+F3T69OlYunQpDh06hJiYGIwdOxb29vYYPHgwAPmPxd/fH++99x7+/vtvnDt3DlOmTMHIkSOFw63Ro0dDKpUiKCgIcXFx2Lt3L7766ivRYbk6dXNuCjszfTzrQEsCwM5MH92cK89XXlPGxsYICgrC7NmzcerUKcTGxmLcuHHP3Wvx9fWFl5cXBg8ejNDQUNy6dQvnz5/H/Pnzhftz58+fj+zsbKxYsQJz5syBi4sL3n33XVE7+vr6CAwMxJUrVxAeHo5p06Zh+PDhwpFAnz59kJeXh7i4uGqt08yZM3Hy5EksWbIEiYmJ2LFjB77++mvMmjWrmlunai1atIBUKsX69etx8+ZNHDp0CEuWLHn+TAYG+Cs2Fnq6usIOANNMGh08u3btioMHD2L37t1o3749lixZgrVr1yIgIECoM2fOHEydOhUTJ05E165dkZeXh2PHjonOlezcuRNubm7o27cvBgwYgJ49e4ru4TQzM0NoaCiSk5PRuXNnzJw5EwsXLhTdC6pO2loSLBokP33xdABVvF806CWVXixSCAkJgbe3NwYNGgRfX1/07NkTnTt3fmZ9iUSCI0eOoFevXhg/fjxcXFwwcuRI3L59GzY2Njhz5gzWrl2LHTt2CLl0fvjhB4SHh2Pjxo1CO61bt8Zbb72FAQMGoH///ujQoQM2bNggTLewsMCQIUOqfTP5yy+/jH379mHPnj1o3749Fi5ciM8//7xWRqS3srLC9u3b8dNPP+Gll17CihUrsHr16ufPZGCA3cePI8DfH4YqykPPagcngFORukgAdyw2BcG/XRXdrmRnpo9Fg16Cf3u758ypeTgBXOUePXwI1zZt8M+OHXDu319+DlRFOAGcahPAafQFIybm394O/V6yxd/JGUjPLYS1ifxQvTb2OOuD8gng3N3d1d0dlbh1+zY2LFwI52bNgIIClQZPplocPOsZbS0JvFqp7nak+q6hJYDr0qULulhZAQ8fAk+eAE1Vex6bqY7m7ruzRqvRJ4BT7G0WFKi3H+y5OHgypmk4eNYLHDwZ0zSK4FlcDKjw8USmWhw8GdM0OjqAYpAa3vvUWBw8GdNEfOiu8Th4MqaJFCNzqXgkIKY6HDwZ00S856nxOHiyaqssAVxCQgJsbW01fqzI51m8eLEwGhcAzJ07F1OnTq1yvu3bt8Pc3PyZ7dSIYs+zoABogA8BPr3N6iMOnkwl5s2bh6lTp8LExETdXVGZWbNmYceOHbh582a15ys/0leN6OkBEgkgkwFFRS/W1lPGjx8vDJzDao6DZ30jKwOSw4GY/fJ/ZWXq7hHu3LmD33//vU6e9lF1KoXnsbS0hJ+fn2jAEmUYGxvDwuIFnwKTSF740P1FE9fV5baujzh41idXDwFr2wM7XgcOBMn/XdteXl5L8vPzMXbsWBgbG8POzg5r1qypUGffvn3w8PBAs2bNhDLFYdnx48fRtm1bGBsbw9/fHykpKUIdmUyGJUuWoHnz5tDT00PHjh1x7NgxYfqtW7cgkUiwd+9e+Pj4QF9fHzt37sS4ceMwePBgLFu2DDY2NjA3N8fnn3+O0tJSzJ49G02bNkXz5s2xbds2UT8/+eQTuLi4wNDQEC1btsSCBQuqDBCDBg3Cnj17qrXNnj5sV/R39erVsLOzg4WFBSZPnixadlFREWbNmoVmzZrByMgInp6eOBMdLZ/45AkeP36MUaNGoVmzZjA0NIS7uzt2794tWm7v3r0xZcoUTJ8+XQj8T1uxYgW+//57/Prrr5BIJJBIJDhz5swztzUAbN26FW3btoW+vj7c3NxEo1sp5vv555/Rp08fYbR9RQ4xhe3bt6NFixYwNDTEkCFD8Pjx42ptU42kkjRyrPazZ8b9SrTIjGiR6VMvM/kr7tcX6P2zffDBB9SiRQs6ceIERUdH0+uvv04mJiai7JlvvPEGvf/++6L5tm3bRrq6uuTr60sXLlygqKgoatu2LY0ePZqI5Jkcv/jiCzI1NaXdu3dTfHw8zZkzh3R1dSkxMZGI/suG6eTkRAcOHKCbN2/SgwcPKDAwkExMTGjy5MkUHx9P//d//0cAyM/Pj7744gtKTEykJUuWkK6uLt29e1fo05IlS+jcuXOUnJxMhw4dIhsbG1q5cqUwXZG1s7xr164RAEpOTn7mNlJk1nxWO4GBgWRqakrvv/8+Xbt2jX777TcyNDSkLVu2CHUmTJhA3bt3p7Nnz9KNGzcoJCSE9PT0KPHAAaLr1+nevXsUEhJCly5doqSkJFq3bh1pa2tTZGSk0IaPjw8ZGxvT7NmzKT4+nuLj40X9LCsro7t379KwYcPI39+fUlJSKCUlhYqKip65rX/88Ueys7MTyg4cOEBNmzal7du3iz4jNzc3+v333ykhIYGGDh1Kjo6OVFJSQkREf/31F2lpadHKlSspISGBvvrqKzI3Nxdts/LqS/ZMDp4qUqvBs6yUaI1bJYGzXABd01ZeT4Vyc3NJKpXSvn37hLLHjx+TgYGBKHh6eHjQ559/Lpp327ZtFdLnfvPNN2RjYyNfpbIysrOzo6VLl4rm69q1K3344YdE9N8Pc+3ataI6gYGB5OjoKEqh6+rqSt7e3sL70tJSMjIyElJQVyYkJIQ6d+4svK8seCo+1zNnzjyzHWWCp6Ojoygl8bBhw2jEiBFEJE81rK2tLaTCVujr40Pzxo0jio6udLkDBw6kmTNnCu99fHyoU6dOz+ynIvXw2LFj6c033xRNe9a2btWqFe3atUtUtmTJEvLy8hLNt3XrVmF6XFwcAaBr164REdGoUaNowIABojZGjBhR74Mnj6pUH9w+D+Q8eE4FAnLuy+s5e6tssUlJSSguLhbSlQDy1BuKVBgKBQUFlY67aGhoKEoEZ2dnJ6ShyMnJQUpKCrp37y6ap0ePHrhy5YqorHzaaIV27dqJxqS0sbFB+/bthffa2tqwsLAQpb3Yu3cv1q1bh6SkJOTl5aG0tLTKcR0N/j3vqMi82K5dO9y+fRsA4O3tjaNHjz53/vL9LZ/mw87ODjExMQCAmJgYlJWVwcXFRTRPUVERLHR1gaIilBUXY9nKldi3bx/u37+P4uJiFBUVVcjU+ryBqpVRflvn5+cjKSkJQUFBeO+994Ty0tJSmJmZiebr0KGDaN0AID09HW5ubrh27RqGDBkiqu/l5SU6RVMfcfCsD/LSqq5TnXoqZmlpWWl++6cTbkkkkmrnHAIAIyMjpdqurEyR9CwiIgIBAQEIDg6Gn58fzMzMsGfPnkrP4ZaXkZEBAEJ2zyNHjgjnKg2qMdbm8/qWl5cHbW1tREVFVcijZPxvoA5ZvhxfrV8vZB01MjLC9OnTK1wUqmxbVUf5+RWpnL/99lvRH1AAFfpZfv0UCfwU69dQcfCsD4xtqq5TnXpKKp8ArkWLFgD+SwDn4+Mj1OvUqROuXr1arbZNTU1hZ2eH8+fPo0+fPkL5uXPn0K1bN9WsQDnnz5+Ho6Mj5s+fL5Qp9iCfJzY2Frq6umjXrh0ACMkHValTp04oKytDenq6kG9eQARkZ+PcuXN488038c477wCQB6bExMQK2WWVIZVKUVZW9V0aNjY2sLe3x82bN0Wpb6qrbdu2oiSCgDyhX33HwbM+cOwOmNoDOSmQZ2l/mkQ+3bF7JdNqrnwCOAsLC1hbW2P+/PkVUjj4+flhwoQJKCsrq1YGyqlTp2LFihVo3bo1OnbsiG3btuHy5cvVzkukjDZt2uDOnTvYs2cPunbtisOHD+PgwYNVzhceHg5vb+9q7WVWl4uLCwICAjB27FisWbMGnTp1wsOHD3Hy5El0sLfHwHbt0KZFC+wPDcX58+fRpEkTfPnll0hLS6tR8HRyckJoaCgSEhJgYWFR4RC8vODgYEybNg1mZmbw9/dHUVER/vnnH2RmZiqdIHHatGno0aMHVq9ejTfffBPHjx+v94fsAN+qVD9oaQP+K/9984wUcP4r5PVUTJkEcK+99hp0dHRw4sSJarU9adIkfPzxx5g5cybc3d1x7NgxHDp0CG3atFHlKgAA3njjDXz88ceYMmUKOnbsiPPnz2PBggVVzrdnzx7R+b7asm3bNowdOxYzZ86Eq6srBg8ejAsXLqBFy5YAgM8mTMDLL78MPz8/9O7dG7a2tjW+0X3ChAlwdXVFly5dYGVlVSEV89N1t27dim3btsHd3R0+Pj7Yvn07nJ2dlV7eK6+8gm+//RZfffUVPDw8EBoais8++6xGfdcknABOReoiARyuHgKOfSK+eGTaTB44X3qj5u2qwDfffINDhw7h+PHjStWvD8nIjh49ipkzZyI6Oho6Omo6SCsoAOLiAC0toFMn+c3zNVQftjnACeBYbXjpDcBtoPyqel6a/BynY/da2eOsrkmTJiErKwu5ubkN5hHN/Px8bNu2TX2BExA/pllcLH/PNAIHz/pGS1ultyOpio6OjuhiTEMwdOhQdXdBvsepry/fAy0o4OCpQTR3350xJsdje2okDp6MaToe21MjcfCsQ3xtjtVI+bE9G4H68jvh4FkHFE9fPOHDLlYTij3PwkJAiZvb6zvFU1PVuWdYHfiCUR3Q1taGubm58Jy1oaGh8AhbYyWTyVBcXIzCwkKNvm1GY2hrywNndvZ/e6LVVB+2uUwmw8OHD2FoaKjeuxyUoNm9a0BsbW0BQDRQRWNGRCgoKICBgUGj/0OilOxs+Z4nEWBsXKMm6ss219LSQosWLTS6jwAHzzojkUhgZ2cHa2trHqEb8lHKz549i169elUYNINVYt8+YNs2YMwYoIa3hNWXbS6VSjV2z7g8Dp51TFtbW+PP5dQFbW1tlJaWQl9fX6N/yBrD2Rm4fRv480/5fZ81wNtctTQ/vDPGAMV4mVeuNMhsmvURB0/G6oO2beUXjTIzgfv31d0bBg6ejNUPenqAm5v8/4rEcEytOHgyVl8oDt05eGoEDp6M1RccPDUKB0/G6gsOnhqFgydj9YWHh/zf+HigqEi9fWEcPBmrN+ztgaZN5Y9pXrum7t40ehw8GasvJBLx/Z5MrTh4Mlaf8HlPjcHBk7H6hIOnxuDgyVh9wsFTY3DwZKw+addOnhQuPR1IS1N3bxo1Dp6M1SeGhkCbNvL/896nWnHwZKy+4UN3jcDBk7H6hoOnRuDgyVh9w/d6agQOnozVN4rgefUqwCld1IaDJ2P1jaMjYGIiD5wJCeruTaPFwZOx+qb8Y5p83lNtOHgyVh8pRlji4Kk2HDwZq494z1PtOHgyVh9x8FQ7Dp6M1Uft28v/vX8fePxYvX1ppDh4MlYfmZgALVvK/897n2rBwZOx+ooP3dWKgydj9RUHT7Xi4MlYfcW3K6kVB0/G6ivFnmdsrDwpHKtTHDwZq69atpSP71lYCNy4oe7eNDocPBmrr7S0AHd3+f/50L3OcfBkrD7j4enUhoMnY/UZX3FXGw6ejNVnHDzVhoMnY/WZInjevg1kZ6u3L40MB0/G6jNzc6BFC/n/Y2LU25dGRqODZ1lZGRYsWABnZ2cYGBigVatWWLJkCYhIqENEWLhwIezs7GBgYABfX19cv35d1E5GRgYCAgJgamoKc3NzBAUFIS8vT1QnOjoa3t7e0NfXh4ODA1atWlUn68jYC+NDd7XQ6OC5cuVKbNy4EV9//TWuXbuGlStXYtWqVVi/fr1QZ9WqVVi3bh02bdqEyMhIGBkZwc/PD4WFhUKdgIAAxMXFISwsDL///jvOnj2LiRMnCtNzcnLQv39/ODo6IioqCiEhIVi8eDG2bNlSp+vLWI1w8FQP0mADBw6kd999V1T21ltvUUBAABERyWQysrW1pZCQEGF6VlYW6enp0e7du4mI6OrVqwSALly4INQ5evQoSSQSun//PhERbdiwgZo0aUJFRUVCnU8++YRcXV2V7mt2djYBoOzs7OqvaCNUXFxMv/zyCxUXF6u7K/Xfnj1EAJGX13OrNfZtrurfqEbveXbv3h0nT55EYmIiAODKlSv4888/8dprrwEAkpOTkZqaCl9fX2EeMzMzeHp6IiIiAgAQEREBc3NzdOnSRajj6+sLLS0tREZGCnV69eoFqVQq1PHz80NCQgIyMzNrfT0ZeyHl9zxlMvX2pRHRUXcHnmfu3LnIycmBm5sbtLW1UVZWhi+++AIBAQEAgNTUVACAjY2NaD4bGxthWmpqKqytrUXTdXR00LRpU1EdZ2fnCm0opjVp0qRC34qKilBUVCS8z8nJAQCUlJSghNPBVkmxjXhbqYCTE3T09CDJz0dJYiLQqlWl1RTbuqiwEPF/HUdBVgoMzO3QpktfaOtodChQCVV/1zR6i+3btw87d+7Erl270K5dO1y+fBnTp0+Hvb09AgMD1dq35cuXIzg4uEJ5aGgoDA0N1dCj+iksLEzdXWgQfJo3h3lSEi7t2IGUV155bt2Tp079+z9rIKMMN0NDa7+DGuDJkycqbU+jg+fs2bMxd+5cjBw5EgDg7u6O27dvY/ny5QgMDIStrS0AIC0tDXZ2dsJ8aWlp6NixIwDA1tYW6enponZLS0uRkZEhzG9ra4u0tDRRHcV7RZ2nzZs3DzNmzBDe5+TkwMHBAf3794epqemLrHajUFJSgrCwMPTr1w+6urrq7k69p33wIJCUhM66upANGFBpncsnduFBsTl8o6dBSv9dUJX9e/NKrOcqdOg7ui66qxaKo0NV0ejg+eTJE2hpiU/LamtrQ/bveR1nZ2fY2tri5MmTQrDMyclBZGQkPvjgAwCAl5cXsrKyEBUVhc6dOwMATp06BZlMBk9PT6HO/PnzUVJSIvyQw8LC4OrqWukhOwDo6elBT0+vQrmuri4Hg2rg7aUiHTsCO3ZAOy4O2pVsz7LSUthfWIEHHisgpULoygpF02UENPsrGFq+AQ32EF7V3zONvmA0aNAgfPHFFzh8+DBu3bqFgwcP4ssvv8SQIUMAABKJBNOnT8fSpUtx6NAhxMTEYOzYsbC3t8fgwYMBAG3btoW/vz/ee+89/P333zh37hymTJmCkSNHwt7eHgAwevRoSKVSBAUFIS4uDnv37sVXX30l2rNkTKNVMUBIfORxWCPjmbNrSQBbPEZ85PHa6F2DpNF/YtavX48FCxbgww8/RHp6Ouzt7TFp0iQsXLhQqDNnzhzk5+dj4sSJyMrKQs+ePXHs2DHo6+sLdXbu3IkpU6agb9++0NLSwttvv41169YJ083MzBAaGorJkyejc+fOsLS0xMKFC0X3gjKm0RRD0928CeTlAcbGoskFmfeVakbZegyQEJV7XIfVWE5ODszMzJCdnc3nPJVQUlKCI0eOYMCAAXzYrir29kBKCvDXX8C/p6QU4s4dhsvJd3HEYwsGXJlY4bBdqNdvF9r1GFgXva1zqv6NavRhO2OsGp5z6O7m6Yd0NH3mrDICUmEBN0+/2updg8PBk7GG4jmPaWrr6CC121wA/11dV1C8T/Fa1GAvFtUGDp6MNRRVZNNU3Ib0SCLeA02XWOBK93Xo5Kfee6frG/4zw1hDUX7PkwiQSCqtZjH7IuIunkJB5n0YNGkGN08/2PIeZ7XxFmOsoXB1BXR15YMi37373zifT9HW0WmwF4XqEh+2M9ZQSKVA27by//PwdLWOgydjDQmP7VlnOHgy1pBwKuI6w8GTsYaE9zzrDAdPxhoSRfBMTAQKCtTblwaOgydjDYmtLWBpKR9R/upVdfemQePgyVhDIpFUebM8Uw0Onow1NHzes05w8GSsoeHgWSc4eDLW0JS/XYlHnKw1HDwZa2heegnQ0gIeP5aP78lqBQdPxhoafX35c+4AH7rXIg6ejDVEfN6z1nHwZKwh4uBZ6zh4MtYQcfCsdRw8GWuIFDfKX7sGFBerty8NFAdPxhqi5s0Bc3OgtFQeQJnKcfBkrCGSSPjQvZZx8GSsoeLgWas4eDLWUHHwrFUcPBlrqDh41ioOnow1VO3ayc99pqYC6enq7k2Dw8GTsYbK2Bho1Ur+/5gY9falAeLgyVhDprjfU4mEcGWlpYg7dxj//L4FcecOo6y0tJY7V7/pqLsDjLFa1KEDcOBAlec9Lx3fAfuIYLTDY6EsLcwCD7wWoZNfYG33sl7i4MlYQ6bERaNLx3fA4/w0+RvJf+VW9BhW56fhEsABtBJ82M5YQ6YInnFx8qeNnlJWWgr7iGAAgJZEPE3x3i4imA/hK8HBk7GGzMlJfuGouBhISKgwOT7yOGzwuELgVNCSALZ4jPjI47Xbz3qIgydjDZmWFuDuDgCQxMZWmFyQeV+pZpSt15hw8GSsofv30F1Sye1KBk2aKdWEsvUaEw6ejDV0zwmebp5+SIMFZM/IEycjIBUWcPP0q80e1kscPBlr6P6917Oy4Kmto4MHXosAoEIAVbxP8VoEbR2+MedpHDwZa+jatwcASO7dg25uboXJnfwCcaX7OjyUWIjK0yUWuNJ9Hd+m9Az854Sxhs7MTH7V/dYtmN6+XWmVTn6BKOsbgLjI4yjIvA+DJs3g5ukHW97jfCbeMow1Bh06yIPnrVvPrKKto4N2PQbWXZ/qOT5sZ6wx+Pei0bP2PFn1cfBkrDH4N3iaPWfPk1UPB0/GGoN/g6fJ7dtAWZmaO9MwcPBkrDFo3RpkYACd4mIgKUndvWkQOHgy1hhoa4PatQNQ+f2erPo4eDLWWCiecefgqRIcPBlrJIiDp0px8GSskaDnjK7Eqo+DJ2ONBCke00xOBnJy1Nyb+o+DJ2ONhYUFCiz+fX6dD91fGAdPxhqRbCcn+X+qSAjHqsbPtjPWiOQ4OcE2Kkqp4FlWWor4pwYK4aHp/sNbgrFGJMfRUf4fTkX8wjh4MtaI5CgO22NiAJlMnuPoKZyKWDl8zpOxRiTP3h4klQK5uUAlIyxxKmLlcfBkrBEhHR2gbVv5m0oO3TkVsfI4eDLWyNC/IyzhypUK0zgVsfI4eDLWyCieNKpsz5NTESuPgydjjczzgienIlYeB0/GGhkheN64AeTni6ZxKmLlcfBkrLGxtgZsbAAiIC6uwmRORawc/vPBWGPUoQMQFiY/dO/WrcJkTkVcNd4SjDVG5YPnM3Aq4ufjw3bGGqPn3K7ElMPBk7HGyMND/m90tPzcJ6s2Dp6MNUZuboCODpCVBdy7p+7e1EscPBlrjPT05AEU4LE9a4iDJ2ONleK8JwfPGlFr8Dx79iwG25CldwAAIABJREFUDRoEe3t7SCQS/PLLL6LpRISFCxfCzs4OBgYG8PX1xfXr10V1MjIyEBAQAFNTU5ibmyMoKAh5eXmiOtHR0fD29oa+vj4cHBywatWqCn356aef4ObmBn19fbi7u+PIkSOqX2HGNAkHzxei1uCZn58PDw8PfPPNN5VOX7VqFdatW4dNmzYhMjISRkZG8PPzQ2FhoVAnICAAcXFxCAsLw++//46zZ89i4sSJwvScnBz0798fjo6OiIqKQkhICBYvXowtW7YIdc6fP49Ro0YhKCgIly5dwuDBgzF48GDEcpZB1pBx8HwxpCEA0MGDB4X3MpmMbG1tKSQkRCjLysoiPT092r17NxERXb16lQDQhQsXhDpHjx4liURC9+/fJyKiDRs2UJMmTaioqEio88knn5Crq6vwfvjw4TRw4EBRfzw9PWnSpElK9z87O5sAUHZ2ttLzNGbFxcX0yy+/UHFxsbq70mhU2Ob37hEBRNraRAUF6u1cHVD1b1Rjb5JPTk5GamoqfH19hTIzMzN4enoiIiICI0eOREREBMzNzdGlSxehjq+vL7S0tBAZGYkhQ4YgIiICvXr1glQqFer4+flh5cqVyMzMRJMmTRAREYEZM2aIlu/n51fhNEJ5RUVFKCoqEt7n/JvKtaSkBCUlJS+8/g2dYhvxtqpdZaWluP7PSRRkpUDPzBZAuW1uZQWdpk0hychASXQ00KmTGnta+1T9XdPY4JmamgoAsLGxEZXb2NgI01JTU2FtbS2arqOjg6ZNm4rqODs7V2hDMa1JkyZITU197nIqs3z5cgQHB1coDw0NhaGhoTKryACEhYWpuwuNhDWQKQMg3ubdmzWDVUYGYn78EXdTUtTVuTrx5MkTlbanscFT082bN0+0t5qTkwMHBwf0798fpqamauxZ/VBSUoKwsDD069cPurq66u5OgxN9chfaR84B8F/6jBItfYS5r4Nv9DTEd/0cHfqOhtbJk0BMDDy0tOA+YIAae1z7FEeHqqKxwdPWVn6IkZaWBjs7O6E8LS0NHTt2FOqkp6eL5istLUVGRoYwv62tLdLS0kR1FO+rqqOYXhk9PT3o6elVKNfV1eVgUA28vVSvrLQUzf8Khi4VygPnUw8Q6cgK0eyvYGj5BkD739+SdmwstBv456Dq75nG3ufp7OwMW1tbnDx5UijLyclBZGQkvLy8AABeXl7IyspCVFSUUOfUqVOQyWTw9PQU6pw9e1Z0viMsLAyurq5o0qSJUKf8chR1FMthrD6pVh6i8s+482Oa1aLW4JmXl4fLly/j8uXLAOQXiS5fvow7d+5AIpFg+vTpWLp0KQ4dOoSYmBiMHTsW9vb2GDx4MACgbdu28Pf3x3vvvYe///4b586dw5QpUzBy5EjY29sDAEaPHg2pVIqgoCDExcVh7969+Oqrr0SH3B999BGOHTuGNWvWID4+HosXL8Y///yDKVOm1P1GYewFVSsPUbt28vTDjx4BTx19sSqo5Jp9DZ0+fZogP6gQvQIDA4lIfrvSggULyMbGhvT09Khv376UkJAgauPx48c0atQoMjY2JlNTUxo/fjzl5uaK6ly5coV69uxJenp61KxZM1qxYkWFvuzbt49cXFxIKpVSu3bt6PDhw9VaF75VqXr4VqXaE/vn70SLTCu8ioOt5ds82Jpokam8HhGRq6v8lqXjx9Xb8Vqm6t+ohIj31VUhJycHZmZmyM7O5gtGSigpKcGRI0cwYMAAPuepYmWlpXi01AVWJD50L9HSxxGPLfC/PBGZZASrzxLl6TSGDwd++glYtQqYPVt9Ha9lqv6N1viw/c6dOwgPD8fx48dx8eJF0T2PjDH1eV4eIgVRHiJ+0qhGqnW1/datW9i4cSP27NmDe/fuofxOq1Qqhbe3NyZOnIi3334bWloaey2KsQavk18gLgGwjwiGDR6LpsV6rkLn8nmIyo/tyZSmdISbNm0aPDw8kJycjKVLl+Lq1avIzs5GcXExUlNTceTIEfTs2RMLFy5Ehw4dcOHChdrsN2OsCp38AmH5WSLi+u3CP11CEN/n/wAAHfqOFldU7HleuwYUF9dxL+svpfc8jYyMcPPmTVhYWFSYZm1tjVdffRWvvvoqFi1ahGPHjuHu3bvo2rWrSjvLGKue8nmISkpKcLOy0cJatABMTf+/vTuPb6pO9wf+SfeyNF1omxbZGVsodAEEyyCbQFkuDDjjAogwog6OO1zkcq+CzHh/DvryqjM6IHoVVJBxo3dYBCKFKUIBoZTSshYKKDQttrSlLG2anN8fpyc2bdMmaZKTc/J5v159kSbfnj6J+HDO9zzf7wNUVwOnTwNSa2Jqld3J87XXXrP7oBMnTnQqGCKSgUYjnn1+/7146d4keZrq63GqSRdN9m138obR2LFjUVlZ2ez56upqjB07tt1BEZGH2bhpdHTHOvz86p1I0s/CkMOLkaSfhZ9fvRNHd6yTIUjv4tQ/H3v27EFdC3Mjt2/fxt69e9sdFBF5WAvdNI/uWIeU/c+K3zQqeYoWyhG9/1kchTiv6qscSp75jf5VOnHihNWuQyaTCdu3b0fXrl1dFx0ReUaTM09TfT3ic8Rdw5ou8/TTiCVQcTkrYLp3ts9ewjv0rlNTU6HRaKDRaFq8PA8NDcXf/vY3lwVHRB4izXOWlABXr+LUmUNIQrnVGWdj0vr4woM7LDekfI3dybO6uhrnz58HAPTu3RuHDh1CdHS05fWgoCDExMTA39/f9VESkXt16gT06QOcOwccP45bNx1YH++j7E6eERERKCkpQUxMDEaNGoW+ffsiPDzcnbERkSclJ4vJMz8foXf9yq4fCY3w3Wk6u++2d+rUCeXl4kqFplu8EZEKNJr3TByWgVJE2VzeaRYAA6KQOCzDc/F5GbvPPMeNG4cxY8agX79+EAQBM2bMsOoL1FhWVpbLAiQiD2mUPKX18dH7n4VZsL5pJCXUkvTl0PnozSLAgeT52WefYd26dTh37hz+9a9/ISkpib16iNRESp6FhUB9vc318WWaKJSkL/fpMiXAgeQZGhqKBQsWAAAOHz6MlStXcs6TSE169wY6dgRu3ACKioDERKRlzIXp3tkobLLCyJfPOCVOfQK7d+92dRxEJDc/P7Fk6cABsd4zMRGA9fp4LtX8hd03jP7yl7/Y3brz4MGD2Lp1q9NBEZFMWtnbk0s1rdmdPE+cOIEePXrgj3/8I7799ltcvXrV8lp9fT3y8/Px97//HcOHD8eDDz6Izp07uyVgInKjVta4p+x/FtGC9d6g0UI5UvY/65MJ1O7k+cknn+C7776D0WjErFmzoNPpEBQUhM6dOyM4OBhpaWn46KOP8Mgjj+DUqVMYOXKkO+MmIndoIXm2tVQTaFiqWV/viQi9hkOTFSkpKfjggw/w/vvvIz8/HxcvXsStW7fQpUsXpKamokuXLu6Kk4g8QVqmefEiUFUFaLU4dXAHl2q2wKmZXj8/P6SmpiI1NdXV8RCRnMLDxc2RL10Szz7vucexVsY+xOnbZGazGUVFRSgrK4PZbLZ6jZfsRAqWnGyVPO1dgulrSzWdSp4HDhzArFmzcPHiRTTtXKzRaGAymVwSHBHJICUF2LLFMu+ZOCwDpfqoZq2MJWZBLJz3taWaTu0kv2DBAgwZMgQFBQWoqKjAtWvXLF8VFRWujpGIPKnJTaPWWhlL318e+pLP1Xs69W7Pnj2Lr776Cn379nV1PEQkNyl5Hj8OmM2An5/NpZrSmegdh17FUX9/n1qy6dSZ57Bhw1BUVOTqWIjIG/TtC4SEiMs0G/bwBcSWGz8NfRmCADSZrfPJek+7zzwbt+B45plnsGjRIhgMBgwcOBCBgYFWY5Olf7mISHkCAoCkJODIEfHSveEK01RfjzsO/RkC2JoDcCB5Si04Gt8gevTRRy2Ppdd4w4hIBZKTf0me990HAKz3bMLu5FlcXOzOOIjIm7TQTZP1ntbsTp49evSwPM7Ozsbw4cMR0OTUvL6+Hvv377caS0QK1PimUQPWe1pz6obRmDFjWixJqqqqwpgxY9odFBHJTEqe584B168DAFtzNOFU8pTmNpsqLy9Hx44d2x0UEcmsSxcI8fEAgJPvv4rCfeIWk23Ve5akL/eJm0WAg3We9zVMHGs0GsybNw/BwcGW10wmE/Lz8zF8+HDXRkhEHnd0xzr073ANwQD67fkrULMapfooIH05jg3/K1tzwMHkqdVqAYhnnp07d0ZoaKjltaCgINx99914/PHHXRshEXmUtHenJsYEFAEoFatnooVyRO9/FseG/xVdXjrj8605HHq3H3/8MQCgZ8+e+Pd//3deohOpTOO9OzWx/uKTZeLGP41rOXHvbJ8oR2qNU3Oey5cvR8eOHVFWVoa9e/di7969KCsrc3VsRORhpw7uQCwaNgCJbUgPpSbLkiKplvPUwR3yBeklnEqe169fx5w5c9C1a1eMGjUKo0aNQteuXfHwww+jqqrK1TESkYdY1Wh28RMzRC2AKsH2OB/lVPJ87LHHcPDgQWzZsgWVlZWorKzEli1bcPjwYfzhD39wdYxE5CFWNZr+GiC60dmnrXE+yqkZ3i1btmDHjh0YMWKE5bmMjAx88MEHmDhxosuCIyLParZ3Z6w/UGoWvxJ8d+/Oljh15hkVFWW5896YVqtFREREu4MiInk027szpiFFlJnsquU01dejcN9WHN6yBoX7tqq6KZxTZ54vvfQSFi5ciE8//RQ6nQ4AYDAYsHjxYrz88ssuDZCIPMtq787YUvFJg7nNWs6jO9YhPmeFuHlIg1J9FK6otP7TqeS5atUqFBUVoXv37ujevTsA4NKlSwgODsbVq1fx/vvvW8bm5ua6JlIi8pi0jLkw3Tsbp7d+joT1j0C4BkQvzIMuLKzF8VJtKACrXZek2tCjDcdUE6eS5/Tp010dBxF5Gf+AACRMexiIXgTN1avwP3MGGDKk2bi2+rqrdZ9Pp97J8uXLXR0HEXkjjUbcJGTXLnF7uhaSp6/u8+nUDSMAqKysxIcffoilS5dadljKzc3F5cus/yJSlRa2p2vMV/f5dOrMMz8/H+PGjYNWq8WFCxfw+OOPIzIyEt988w0uXbqETz75xNVxEpFcmnTTbMpX9/l06sxz4cKFmDdvHs6ePYuQkBDL85MnT0Z2drbLgiMiL9A4eTbt/Abf3efTqeT5ww8/tLiSqGvXrjAYDO0Oiojax6X1lv37A/7+QHk5UFLS7GV7+rqrcZ9Pp95NcHAwqqurmz1/5swZREdHtzsoInKerXrLn+5eDiDc8QOGhAAJCcCJE+LZZ8MmyY3Z6uuu5n0+nTrznDZtGv70pz/BaDQCEDdHvnTpEpYsWYLf/va3Lg2QiOwn1VtGC+VWz0cL5Rhw8EXnD9zGvCcgJtAuL51B4fgNODzkDRSO34Dol86oMnECTibPN998EzU1NYiOjsatW7cwatQo9O3bF507d8Z///d/uzpGIrJDW/WWjcc5zI7kCYiX8Em/noIh//YEkn49RXWX6o059c60Wi30ej327duHY8eOoaamBoMGDcK4ceNcHR8R2cmeeksAOHt4Fwbe82+OHbxJK2JTfT1ONdlJXs2JsiUOv1uz2Yy1a9fim2++wYULF6DRaNCrVy/odDqbjeGIyP3srresbH7Tp00pKeKfp04hb/OHiDvy/3xmDbstDl22C4KAadOm4bHHHsPly5cxcOBAJCUl4eLFi5g3bx5mzJjhrjiJqA1211uGxzl+8K5dgYgIoL4eyVueb3FONWX/szi6Y53jx1Yoh5Ln2rVrkZ2djV27duHo0aP4/PPPsXHjRhw7dgzfffcdsrKyWCBPJBN76i0B4FdD7nX84BoNhAEDxMdlJptzqnE5K1S9DV1jDiXPzz//HP/5n/+JMWPGNHtt7Nix+I//+A+sX7/eZcERkf3sqbeUxjmjIkbcw9ev1Nzi677W38ih5Jmfn9/qTvGTJk3CsYYJZSLyvLSMuTg2/K+4qomyer5ME4WCYa+369jX4xpqRJu05GhKbWvYbXHon6CKigrExsbafD02NhbXrl1rd1BE5DxpL86mfdWjBAE/bdvm9HHNyWkAPhNbcrRCbWvYbXEoeZpMJgS0csrv7++Peh+Z7yDyZlK9ZWPmhkUtzurxwHwITyyC5oYA1JiBTtYXrr7W38ih5CkIAubNm4fg4OAWX6+trXVJUETkOfbWbPprtbjdNRYhl0thLjXDr1HybLyGXecj9Z4Ovcu5c9uu4XrkkUecDoaIPMvRvkMhd48Avv4aN8pC0LnPL1eZal7DbotDyfPjjz92VxxE5GFO9R1KTga+/hodY6egcPyDVmervnLGKfGtd0tEANrRd6hhmaZfQQGSfu3bNd1Ot+EgIuU6dXAHYlHeLHFKbNZsSmvcCwuBdt6AUjomTyIf5HTfoZ49gU6dgLo64OxZ1wemIEyeRD7I6b5Dfn7NdljyVUyeRD6oXX2H7NzbU+1kTZ7Z2dmYOnUq4uPjodFokJmZaXnNaDRiyZIlGDhwIDp27Ij4+Hg88sgjuHLlitUxKioqMHv2bISFhSE8PBzz589HTU2N1Zj8/Hzcc889CAkJQbdu3fD6682XqX355ZdITExESEgIBg4ciG3tWIlB5O3a1XeIyROAzMnzxo0bSElJwXvvvdfstZs3byI3Nxcvv/wycnNz8c033+D06dOYNm2a1bjZs2ejsLAQer0eW7ZsQXZ2Np544gnL69XV1ZgwYQJ69OiBI0eO4I033sArr7yCNWvWWMbs378fM2fOxPz583H06FFMnz4d06dPR0FBgfvePJHMWlsHf2z4X23XbDJ5igQvAUDYtGlTq2MOHTokABAuXrwoCIIgnDhxQgAg/PDDD5Yx3377raDRaITLly8LgiAIf//734WIiAihtrbWMmbJkiVCQkKC5fsHHnhAmDJlitXvGjZsmPCHP/zB7virqqoEAEJVVZXdP+PL6urqhMzMTKGurk7uUHyGrc+83mgUCr7fIvyw+X2h4PstQr3R2PqBKisFQWxCLAjl5W6M2LVc/f+oouo8q6qqoNFoEB4u7u6Sk5OD8PBwDBkyxDJm3Lhx8PPzw8GDBzFjxgzk5ORg5MiRCAoKsozJyMjAypUrce3aNURERCAnJwcLFy60+l0ZGRlW0whN1dbWWi1HlbqJGo1GS2M8sk36jPhZeU5rn/mdQydYHpsFwWodvKm+HmcP78KtyhKEhsfhV0PuRXDPntBcuID6o0chjBzp/uBdwNV/1xSTPG/fvo0lS5Zg5syZCAsLAwAYDAbExMRYjQsICEBkZKSlf7zBYECvXr2sxkg7QxkMBkRERMBgMDTbLSo2NrbVHvSvvfYaVqxY0ez5nTt3okOHDo6/QR+l1+vlDsHnOP+ZxwAVJpzfuRNDY2IQd+ECTm7ciPNN7jF4q5s3b7r0eIpInkajEQ888AAEQcCqVavkDgcAsHTpUquz1erqanTr1g0TJkywJHeyzWg0Qq/XY/z48QgMDJQ7HJ/g6Geev2uDpV1x42J6swD4h4gJM8lkQuLkyW6J19Wkq0NX8frkKSXOixcvIisryyox6XQ6lJWVWY2vr69HRUUFdDqdZUxpaanVGOn7tsZIr7ckODi4xd2lAgMDmQwcwM/L8+z5zE319bjjwAoECrfFxNn0jnyM+ITm+HHF/PdzdZxeXecpJc6zZ8/iu+++Q1SU9V3B9PR0VFZW4siRI5bnsrKyYDabMWzYMMuY7Oxsq/kOvV6PhIQEREREWMbs2rXL6th6vR7p6enuemtEXq3N5Zs6MXUIx48DptZ3ljfV16Nw31Yc3rIGhfu2qqbHkaxnnjU1NSgqKrJ8X1xcjLy8PERGRiIuLg6/+93vkJubiy1btsBkMlnmICMjIxEUFIR+/fph4sSJePzxx7F69WoYjUY8/fTTeOihhxAfHw8AmDVrFlasWIH58+djyZIlKCgowDvvvIO33nrL8nufe+45jBo1Cm+++SamTJmCjRs34vDhw1blTES+pM3lmxF+QADgV1sLFBUBCQktDnN0yzslkfXM8/Dhw0hLS0NaWhoAYOHChUhLS8OyZctw+fJl/POf/8RPP/2E1NRUxMXFWb72799vOcb69euRmJiIe++9F5MnT8aIESOskp5Wq8XOnTtRXFyMwYMHY9GiRVi2bJlVLejw4cOxYcMGrFmzBikpKfjqq6+QmZmJAVK3QCIf0+byTT8NENOQPo4fb3GItOWdWtsUy3rmOXr0aAiCjfVhQKuvSSIjI7Fhw4ZWxyQnJ2Pv3r2tjrn//vtx//33t/n7iHxB4rAMlOqjEC20fOluFoDbsR3Q4UqNuMb9d7+zet3pLe8UxKvnPIlInjlDe5Zvlt/9G/FBCxuEOL3lnYIoM+UT+Qg55wzTMubiKID4nBWIbfT7LS03QnoCq9a3eNnu9JZ3CsLkSeSlnGqT4WK22hjrAgKAigpx0IULQFUVoNVafs7pLe8UhJftRF6orTlDoGHO0EOX8Em/noIh//YEkn495Zc5yshI4I47xMdNNtFp15Z3CsHkSeSFFDNnaGNj5HZteacQTJ5EXkgxc4ZS8mxh3tPpLe8UQrlpn0jFFDNn2Mbenq3OmSqc8t8BkQrZU2dZpvGCOcPGZ55ms9jjqAlpzlRteNlO5IUUM2d4551AUBBw/bp4192HMHkSeSlFzBkGBgL9+4uPfawtBy/bibyYIuYMU1KAvDwxeU6fLnc0HuNF/wWIqCVeP2foow3heNlORO3D5ElE5AQpeRYVATduyBuLBzF5ElH7xMQAsbFiM+ImyzTVjMmTiNrPBy/dmTyJnGQyC8g5V47/y7uMnHPlMNnaBcMXpKSIf/pQ8uTddiInbC8owYrNJ1BSddvyXJw2BMun9sfEAXEyRiYTnnkSUVu2F5Tgyc9yrRInABiqbuPJz3KxvaBEpshk1Dh52tE+Rw2YPIkcYDILWLH5RNM25gB+aW2+YvMJ37uET0wEAgKAykrgp5/kjsYjmDyJHHCouKLZGWdjAoCSqts4VFzhuaC8QXCwmECBFnsaqRGTJ5EDyq7bTpzOjFMVH5v3ZPIkckBM5xCXjlMVJk8ismVor0jEaUNgozsGNBDvug/tFenJsLyDj5UrMXmST3K2RtPfT4PlU8Ut2JomUOn75VP7w99W8yE1k848T58Gbqt/2oJ1nuRz2lujOXFAHFY9PKjZMXS+XOcJAHFxQFQUUF4OnDgBDBokd0RuxeRJPkWq0Wx6ninVaK56eJDdCXR8fx0OFVeg7PptxHQWL9V98oxTotGIZ5+7d4uX7kyeRMpiMgstJrW2ajQ1EGs0x/fX2ZUE/f00SO8T1eY4n9I4eaockyepSmuX5NrQILtrNJkUnWSjj7sa8YYRqUZbyyb1Jwx2HccnazRdxYeWaTJ5kirYs2zy//Ku2HUsn6zRdJWkJLH98M8/Awb7/rFSKiZPUgV7lk2W36hDZMdA1mi6U2io2I4YUP28J5MnqYK9l9ozUrsCYI2mW/nIvCeTJ6mCvZfa4/rrsOrhQdBprcfrtCF2lylRG6Tkefy4vHG4Ge+2kypIyyYNVbdbnPfUQEyQUtkSazTdyEfWuDN5kipIyyaf/CwXGsAqgbZ0Sd5WjaatWlGyg5Q8T54E6uqAoCB543ETJk9SDVctm2SLjXbq3h3QaoGqKuDUqV+SqcoweZKqtHfZpKuWb/o0aZnm3r3ipTuTJ5EyOLts0tXLN31a4+SpUrzbTtSALTZcyAduGjF5EjVgiw0X8oFaTyZPogZsseFCAwaIfxoMwNWr8sbiJkyeRA3YYsOFOnUC+vQRH6u0WJ7Jk1TD2dYaErbYcDGVz3vybjupgqtqM9liw4WSk4FNm5g8ibyVq2sz2WLDRVR+04jJkxTNXbWZbLHhAlLyLCwE6uuBAHWlG855kqKxNtOL9e4NdOwI1NYCZ87IHY3LMXmSorE204v5+QEDB4qPVXjHncmTFI21mV5OxXfcmTxJ0Qb3iEBbU5l+GnEcyYDJk0hetmo4j1y8hrbKOc2COI5koOI77uq6/UWq1FoNZ2292a5jcM5TJlLy/PFH4No1IEI9VwA88ySv1lYv9gs/37TrOJzzlIlWC/ToIT5W2U0jJk/yWvb0Yt/4wyXowoK5Ht2bSXfcVTbvyeRJHif1BwLEOk1ba9DtreGcObQ7AK5H91oqvWnE5Eketb2gBCNWZuHRdT8AAB5d9wNGrMzC9oKSZmPtnafs2aUj2wl7M5UmT94wIo9pvAY92P+X522tQXekhjO9TxTXo3urlBTxz+PHAbNZLJ5XASZP8ghn1qA70osd4Hp0r9W3LxASAty8CZw/L36vAur4J4C8msksYO2+YofXoHN/TZUICACSksTHKqr3ZPIkt5LmOP+89aRd45vOc0r7a3I+U+FUOO/Jy3ZyG1v7bLampXlO7q+pAtK8J5MnUetam+NsSdP5y6Y4n6lwKjzz5GU7uUVbNZot4fylikmF8ufPA9evyxuLizB5kls4spbcTwM8MbIX5y/VrEsXID5efFxQIG8sLsLkSW7hyFpyQQDWZBe3WChPKqKyS3cmT3KLtnqgNybNi67YfMLhdsGkIEyerpOdnY2pU6ciPj4eGo0GmZmZNscuWLAAGo0Gb7/9ttXzFRUVmD17NsLCwhAeHo758+ejpqbGakx+fj7uuecehISEoFu3bnj99debHf/LL79EYmIiQkJCMHDgQGzbts01b1Jl7O2N3lqNZkvYa8gHqGxvT1mT540bN5CSkoL33nuv1XGbNm3CgQMHEC/NmTQye/ZsFBYWQq/XY8uWLcjOzsYTTzxheb26uhoTJkxAjx49cOTIEbzxxht45ZVXsGbNGsuY/fv3Y+bMmZg/fz6OHj2K6dOnY/r06ShQydyMq0g1mzM/OIDnNuZh5gcHbK5LB2zXaLaG+26qmJQ8jx8X52oUTtbFHG44AAAU8UlEQVTkOWnSJLz66quYMWOGzTGXL1/GM888g/Xr1yMwMNDqtZMnT2L79u348MMPMWzYMIwYMQJ/+9vfsHHjRly5cgUAsH79etTV1eGjjz5CUlISHnroITz77LP4n//5H8tx3nnnHUycOBGLFy9Gv3798Oc//xmDBg3Cu+++6543rkBt7avZWgL9fslYvDyln12/h/tuqlhCAhAYCFRXAxcvyh1Nu3n1nKfZbMacOXOwePFiJEnLuxrJyclBeHg4hgwZYnlu3Lhx8PPzw8GDBy1jRo4ciaCgIMuYjIwMnD59GteuXbOMGTdunNWxMzIykJOT4463pTj27KvZ2nylv58G837dq9U5UO676QOCgoD+4lSOGuY9vbpIfuXKlQgICMCzzz7b4usGgwExMTFWzwUEBCAyMhIGg8EyplevXlZjYmNjLa9FRETAYDBYnms8RjpGS2pra1FbW2v5vrq6GgBgNBphNBrtfIfKcKi4AhU1t6x2QmqqouYWDhSVtZr8lk1JwAv/yAMABPmJiTbYT7Ak1GVTEmA21cNsclXk1Jj091LOv5/+AwbA79gxmHJzYZ40yaO/29Xv22uT55EjR/DOO+8gNzcXGo33FU6/9tprWLFiRbPnd+7ciQ4dOsgQkXu9PrTtMT+fPIBtbSxhX9nkOH8e8ksPorriI9hW7ERw5BC9Xi/b7+4TGIgBAAx6PQ6npXn0d9+8aV/LFnt5bfLcu3cvysrK0L17d8tzJpMJixYtwttvv40LFy5Ap9OhrKzM6ufq6+tRUVEBnU4HANDpdCgtLbUaI33f1hjp9ZYsXboUCxcutHxfXV2Nbt26YcKECQgLC3PiHXuvQ8UVls2LW/PR3Lvsuuw2mQX8cP4qKs4cRuSdQ3BX72iuLPIAo9EIvV6P8ePHN7t/4CmawEBg7VrE//wzJk+e7NHfLV0duorXJs85c+a0OA85Z84c/P73vwcApKeno7KyEkeOHMHgwYMBAFlZWTCbzRg2bJhlzH/913/BaDRa/sLo9XokJCQgoqGTX3p6Onbt2oXnn3/e8rv0ej3S09NtxhccHIzg4OBmzwcGBsr2F9Nd7u4bg8hOoTb31QSA8NBAwM8ffv4BbSbCwIZjbjsj/qm2z8vbyfp3dNAgAICmqAiBRiPgwas0V79nWW8Y1dTUIC8vD3l54jxYcXEx8vLycOnSJURFRWHAgAFWX4GBgdDpdEhISAAA9OvXDxMnTsTjjz+OQ4cOYd++fXj66afx0EMPWcqaZs2ahaCgIMyfPx+FhYX4xz/+gXfeecfqrPG5557D9u3b8eabb+LUqVN45ZVXcPjwYTz99NOe/1C8kD01m5W3jJj94cFWS5eIEBsLREeLO8oXFsodTbvImjwPHz6MtLQ0pDXMfSxcuBBpaWlYtmyZ3cdYv349EhMTce+992Ly5MkYMWKEVQ2nVqvFzp07UVxcjMGDB2PRokVYtmyZVS3o8OHDsWHDBqxZswYpKSn46quvkJmZiQEDBrjuzSqcvTWbbZUukY/TaFSzPZ2sl+2jR4+G4ECx7IULF5o9FxkZiQ0bNrT6c8nJydi7d2+rY+6//37cf//9dsfii6R9NQ+cK8dTG3JReav53UtbLTWILJKTge++U3zy9Oo6T/I+/n4a+PlpWkycEi61pFapZI07kyc5zN4llFxqSS1qvMZdwcs0mTzJYY60BCZqpl8/wN8fuHYNaFhGrURMnuSwtrab41JLalVIiLjOHVD0pTuTJ1m4Yrs5tgQmu6hg3tNri+TJs7YXlGDF5hNWuybFaUOwfGr/FttjSKVLTX9G18rPEFkkJwMbNzJ5krLZahEs1Wza6o/OlsDkNBVsjMzk6ePa2m6urZpNtgQmp0jJ89QpoLYWaGGps7fjnKePa6tFMGs2yS3uuAOIiABMJuBkG1txeSkmTx/Hmk2ShUaj+Et3Jk8fx5pNkk3jnkYKxOTp41izSbJReLkSk6fKtVW7yZpNko3CL9t5t13F7K3dZM0mySIpSZz7LCsDSkvFvT4VhMlTpRyt3WTNJnlcx47Ar34FnDkjXrqPHy93RA7hZbsKOdsqWKrZ/E1qV6T3iWLiJPdT8Lwnk6cKsXaTFEPB855MnirE2k1SDJ55kjdh7SYphpQ8T5wAjLa7E3gjJk8VYu0mKUbPnkDnzmLiPH1a7mgcwuSpQs7Wbtq7nyeRyzRepqmwS3eWKqmUo7Wbju7nSeQyycnAvn1i8pw1S+5o7MbkqWL21m46u58nkUvwzJO8UVv7bbZ3P0+idlNo8uScp8q1NY/JmlCS3YAB4p+XLwPl5fLG4gCeeaqYPfOYrAkl2YWFAb16AcXF4vZ0o0fLHZFdeOapUtI8ZtOzSmkec3tBCQDWhJKXUOClO5OnCjmytp01oeQVmDzJGzgyj2lvTSgA1oCS+yhwjTvnPFXI0XnMtmpCAWDEyizWgJL7SMmzsFBsCufvL288dmDyVCFn5jFt1YTqTxhYA0ru16cPEBoK3LoFFBUBCQlyR9QmXrarkLPzmE338wTg1L6gRA7z9wcGDhQfK2Tek8lThVzVl4g1oORRCpv3ZPJUKWkeU6e1voTXaUPsvtRmDSh5lMJaEXPOU8Xa25eINaDkUQorV2LyVLm21ra3Rpo7NVTdbnHeUwPxTJY1oOQS0pznhQtAVRWg1coaTlt42e5DHN2vkz3dyaMiI4E77hAfK+DSnWeePsLZ/TrZ0508KjkZ+Okn8dJ9xAi5o2kVk6cPaO9+nezpTh6TkgJs26aIeU8mT5Vz1X6d7Zk7JbKbgm4acc5T5Zyp1WQvI5JN43Ils1neWNrAM0+Vc7RWk72MSFZ33gkEBQE1NeJd99695Y7IJp55qpwjtZr27gFK5DYBAUBSkvjYyy/dmTxVzt517oN7RDi0jp2X9uQ2Clmmyct2lZNqNZ/8LBcawCo5Nq7VPHLxmt1zo1W36nhpT+6TkiL+6eXJk2eePsCede72zo1KW9Tx0p7cRiFr3Hnm6SPaqtW0d240M+8K2xSTe0nJ89w58cZRp07yxmMDzzx9SNP9OhsnOHvmRiM7BqLiRp3N43OLOnKJ6GhApwMEQdxZ3ksxeRKAX+ZGbd32EQDMSO1q17G4RR21mwLmPZk8yW5hoUF2jeMWddRuCpj3ZPIkAL8s47RFA2DjD5egCwtmm2JyPwUs02TyJAD2L+OcObQ7AG5RR27WOHkK3llDzORJAOyfp+zZpWO723sQtSkxUVxtVFkJ/Pij3NG0iKVKBMCxZZzpfaK4RR25V1AQ0K+fOOeZnw907y53RM3wzJMAON6uuLWyJyKX8PJ5TyZPlbN3DTpbbpDX8fLkyct2FXN0ezm23CCvwuRJcnC29QZbbpDXkJLn6dPA7dtAiHfVDzN5qlB7W2+w5QZ5hbg4ICoKKC8HTpwABg2SOyIrnPNUIWdabxB5HY3Gqy/dmTxVyNHWG0Rey4s3RmbyVCFHajaJvBrPPMmTHK3ZJPJajXdX8rJlmkyeKsSaTVKN/v0BPz/xppHBIHc0Vpg8Vcqe1htEXi80VGxHDHjdpTtLlVSMNZukCsnJwKlTYvLMyJA7GgsmT5VjzSYpXnIy8MUXXnfmyct2IvJu0k0jJs9fZGdnY+rUqYiPj4dGo0FmZmazMSdPnsS0adOg1WrRsWNH3HXXXbh06ZLl9du3b+Opp55CVFQUOnXqhN/+9rcoLS21OsalS5cwZcoUdOjQATExMVi8eDHq6+utxuzZsweDBg1CcHAw+vbti7Vr17rlPRORg6RypZMngTrbDQg9TdbkeePGDaSkpOC9995r8fVz585hxIgRSExMxJ49e5Cfn4+XX34ZIY3WuL7wwgvYvHkzvvzyS/zrX//ClStXcN9991leN5lMmDJlCurq6rB//36sW7cOa9euxbJlyyxjiouLMWXKFIwZMwZ5eXl4/vnn8dhjj2HHjh3ue/NEZJ9u3QCtFjAaxblPbyF4CQDCpk2brJ578MEHhYcfftjmz1RWVgqBgYHCl19+aXnu5MmTAgAhJydHEARB2LZtm+Dn5ycYDAbLmFWrVglhYWFCbW2tIAiC8OKLLwpJSUnNfndGRobd8VdVVQkAhKqqKrt/xpfV1dUJmZmZQl1dndyh+AxFf+b33CMIgCB8+qnTh3D1/6Nee8PIbDZj69atePHFF5GRkYGjR4+iV69eWLp0KaZPnw4AOHLkCIxGI8aNG2f5ucTERHTv3h05OTm4++67kZOTg4EDByI2NtYyJiMjA08++SQKCwuRlpaGnJwcq2NIY55//nmb8dXW1qK2ttbyfXV1NQDAaDTCaDS65DNQM+kz4mflOUr+zP0GDID/3r0wHT0K84MPOnUMV79vr02eZWVlqKmpwV/+8he8+uqrWLlyJbZv34777rsPu3fvxqhRo2AwGBAUFITw8HCrn42NjYWhoaDWYDBYJU7pdem11sZUV1fj1q1bCA0NbRbfa6+9hhUrVjR7fufOnejQoYPzb9zH6PV6uUPwOUr8zHsASAXw8+7dOLBtm1PHuHnzpktj8trkaTabAQC/+c1v8MILLwAAUlNTsX//fqxevRqjRo2SMzwsXboUCxcutHxfXV2Nbt26YcKECQgLC5MxMmUwGo3Q6/UYP348AgMD5Q7HJyj5M9d06QKsWoUYgwGTJ0926hjS1aGreG3y7NKlCwICAtC/f3+r5/v164fvv/8eAKDT6VBXV4fKykqrs8/S0lLodDrLmEOHDlkdQ7ob33hM0zv0paWlCAsLa/GsEwCCg4MRHBzc7PnAwEDF/cWUEz8vz1PkZ56SAmg00JSUILCyEoiOdvgQrn7PXlvnGRQUhLvuugunT5+2ev7MmTPo0aMHAGDw4MEIDAzErl27LK+fPn0aly5dQnp6OgAgPT0dx48fR1lZmWWMXq9HWFiYJTGnp6dbHUMaIx1DyeztYUTk1Tp1Avr0ER8fPy5vLA1kPfOsqalBUVGR5fvi4mLk5eUhMjIS3bt3x+LFi/Hggw9i5MiRGDNmDLZv347Nmzdjz549AACtVov58+dj4cKFiIyMRFhYGJ555hmkp6fj7rvvBgBMmDAB/fv3x5w5c/D666/DYDDgpZdewlNPPWU5c1ywYAHeffddvPjii3j00UeRlZWFL774Alu3bvX4Z+JKjvYwIvJqyclAUZG4w9LYsXJHI2+p0u7duwWIG5tbfc2dO9cy5n//93+Fvn37CiEhIUJKSoqQmZlpdYxbt24Jf/zjH4WIiAihQ4cOwowZM4SSkhKrMRcuXBAmTZokhIaGCl26dBEWLVokGI3GZrGkpqYKQUFBQu/evYWPP/7YoffibaVK3x6/IvRcskXo0eSrZ8PXt8evyBqfostmFErxn/krr4jlSvPmOfXjrv5/VCMIXrZJnkJVV1dDq9WiqqpK9htGJrOAESuzbLbi0EDcXen7JWNl2yTEaDRi27ZtmDx5svLm3xRK8Z/5pk3AffeJvYyOHHH4x139/6jXznmS89jDiFRJWuNeWAg0WV4tByZPFWIPI1Klnj3FG0e1tcDZs3JHw+SpRuxhRKrk5wcMHCg+9oIdlpg8VYg9jEi1vKghHJOnCrGHEamWF7UiZvJUKfYwIlXyojNPr12eSe3HHkakOtKc548/AteuARERsoXC5Kly7GFEqqLVAj16ABcviss0R46ULRRethORsnhJTyMmTyJSFi+Z92TyJCJlYfIkInKClDyPHwdMJtnCYPIkImXp2xcICQFu3gTOn5ctDCZPIlIWf39gwADxsYyX7kyeRKQ8XjDvyeRJRMrD5ElE5AQvqPVk8iQi5ZGWaZ4/D7i4pbC9mDyJSHmiooCuXcXHBQWyhMDk6WFsBUzkIjLPe3JjEA9iK2AiF0pOBr79VrbkyTNPD9leUIInP8tt1pjNUHUbT36Wi+0FJTJFRqRQMp95Mnl6gMksYMXmE2jpAl16bsXmE7yEJ3JE4+QpQwd1Xra7iNDwH6+6hTt/h85X4HJZ621+L5fdxO78ixja2zf6ChmNRty8eRPV1dXK7CGuQKr7zOPigMBA4Pp1cZ17z56tDpf+3xRclGg1gquO5ON++ukndOvWTe4wiKgNP/74I+644452H4fJ00XMZjOuXLmCzp07Q6Nhm4u2VFdXo1u3bvjxxx8RFhYmdzg+wdc/c0EQcP36dcTHx8PPr/0zlrxsdxE/Pz+X/Gvma8LCwnzyf2Q5+fJnrtVqXXYs3jAiInICkycRkROYPEkWwcHBWL58OYKDg+UOxWfwM3ct3jAiInICzzyJiJzA5ElE5AQmTyIiJzB5EhE5gcmTXC47OxtTp05FfHw8NBoNMjMz2/yZPXv2YNCgQQgODkbfvn2xdu1a9weqIo5+5nv27IFGo2n2ZTAYPBSx8jF5ksvduHEDKSkpeO+99+waX1xcjClTpmDMmDHIy8vD888/j8ceeww7duxwc6Tq4ehnLjl9+jRKSkosXzExMW6KUH24PJNcbtKkSZg0aZLd41evXo1evXrhzTffBAD069cP33//Pd566y1kZGS4K0xVcfQzl8TExCA8PNwNEakfzzxJdjk5ORg3bpzVcxkZGcjJyZEpIt+RmpqKuLg4jB8/Hvv27ZM7HEVh8iTZGQwGxMbGWj0XGxuL6upq3Lp1S6ao1C0uLg6rV6/G119/ja+//hrdunXD6NGjkZubK3doisHLdiIflJCQgISEBMv3w4cPx7lz5/DWW2/h008/lTEy5eCZJ8lOp9OhtLTU6rnS0lKEhYUhNDRUpqh8z9ChQ1FUVCR3GIrB5EmyS09Px65du6ye0+v1SE9Plyki35SXl4e4OHZxtRcv28nlampqrM5giouLkZeXh8jISHTv3h1Lly7F5cuX8cknnwAAFixYgHfffRcvvvgiHn30UWRlZeGLL77A1q1b5XoLiuPoZ/7222+jV69eSEpKwu3bt/Hhhx8iKysLO3fulOstKI9A5GK7d+8WIDYGtfqaO3euIAiCMHfuXGHUqFHNfiY1NVUICgoSevfuLXz88ccej1vJHP3MV65cKfTp00cICQkRIiMjhdGjRwtZWVnyBK9Q3JKOiMgJnPMkInICkycRkROYPImInMDkSUTkBCZPIiInMHkSETmByZOIyAlMnkRETmDyJCJyApMnEZETmDyJiJzA5ElE5IT/D0fVnRfQ/O5cAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 300x600 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "px-NIBrfCuy9",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "\t# Legend Loc -best\n",
        "\t# upper right\n",
        "\t# upper left\n",
        "\t# lower left\n",
        "\t# lower right\n",
        "\t# right\n",
        "\t# center left\n",
        "\t# center right\n",
        "\t# lower center\n",
        "\t# upper center\n",
        "\t# center"
      ],
      "execution_count": 63,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "kjrOGWe6IPiR",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        ""
      ],
      "execution_count": null,
      "outputs": []
    }
  ]
}